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This  report  sunnarizes  the  work  perfonned  by  SRICO  under  contract 
number  DAAL02-90-C-0078  sponsored  by  the  U.S.  Army  Harry  Diamond 
Laboratories.  This  was  a  Phase  I  research  program  entitled 
"Optical  Interferometers  for  Sensing  Electromagnetic  Fields".  The 
overall  goal  of  the  program  was  to  determine  the  feasibility  of 
manufacturing  a  nonmetallic  electromagnetic  field  sensor  to 
measure  wideband  transient  electromagnetic  responses  using 
integrated  optical  waveguide  interferometers.  This  noimetallic 
design  minimizes  the  perturbation  on  the  electric  fields  and  allows 
the  accurate  measurement  of  large  fields. 

The  integrated  optic  device  makes  use  of  a  Mach-Zehnder 
interferometer  in  a  lithium  niobate  electro-optic  crystal  material. 
The  electric  field  sensor  uses  the  electro-optic  properties  of 
lithium  niobate  to  modulate  the  phase  of  the  light  propagating  in 
each  arm  of  the  Mach-Zehnder  interferometer.  This  phase  modulated 
optical  signal  is  converted  to  intensity  modulation  at  the  output 
of  the  interferometer. 

Conventional «  integrated  optic  electric  field  sensors  utilize  a 
dipole  antenna  connected  to  the  metallic  electrodes  on  the 
waveguide  channels  to  apply  voltage  of  opposite  polarities  on  the 
two  channels.  This  is  an  indirect  method  of  measuring  the  electric 


The  need  for  metal  electrodes  could  be  eliminated  if  one  leg  of  the 
interferometer  were  to  be  oppositely  poled.  SRICO  successfully 
demonstrated  a  selective  reverse  poling  of  one  leg  of  the 
interferometer  using  high  temperature  diffusion  of  titanium  into 
lithium  niobate.  A  series  of  experiments  were  performed  to 
detezmine  the  optimum  diffusion  conditions  to  obtain  100%  reverse 
poling  of  one  waveguide  channel  without  affecting  the  rest  of  the 
crystal  substrate.  The  optical  waveguide  is  then  fabricated  by  a 
low  temperature  proton-exchange  process. 

Analysis  of  the  sensor  system  indicates  that  it  can  detect 
electromagnetic  fields  with  a  dynamic  range  greater  than  40  dB  and 
a  signal  bandwidth  of  3  GHz .  The  sensor  could  be  designed  to 
operate  from  DC  to  3  GHz  or  as  a  passband  device  operating  at  a 
frequency  as  high  as  30  GHz. 

The  sensor  is  packaged  with  optical  fibers  for  input  and  output. 
The  optical  fibers  enable  remote  measurement  of  the  field  from  a 
safe  distance  of  100  meters. 

The  electric  field  sensor  developed  during  this  Phase  I  program  has 
been  shown  to  be  feasible  and  a  prototype  sensor  instrument  has 
been  proposed  for  a  Phase  II  SBIR  program. 
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Fiber  optic  technology  has  Bade  great  strides  in  the  last  two 
decades,  especially  for  applications  in  telecomunications  and 
sensors.  Fiber  optics  provides  the  advantages  of  secure 
transnission ,  iiniunity  to  EMI/EMP,  wide  bandwidth,  low  weight  and 
long  link  length  on  the  order  of  several  kiloneters.  Optical 
fibers  can  be  conbined  with  integrated  optical  devices  to  yield 
ainiature,  conpact  and  rugged  conponents  for  sensing  applications. 

In  addition  to  their  conpatibility  with  optical  fibers,  other 
advantages  of  optical  waveguide  devices  include  peraanent  alignaent 
of  coaponents,  coapatibility  with  existing  planar  processing 
technology,  diffraction**free  beaa  propagation,  low  voltage  and/or 
power  requireaents ,  high  aodulation  rates  and  freedoa  froa 
electroaagnetic  interference.  Optical  fibers  can  be  directly 
attached  to  optical  waveguides  without  the  need  for  any  bulk  optic 
coaponents  such  as  lenses  and  polarizers.  Integrated  optical 
sensors  could,  therefore,  be  aade  coapact,  rugged  and  reliable  to 
withstand  severe  ailitary  envlronaents. 

One  popular  integrated  optical  waveguide  device  is  the  Mach-Zehnder 
interferoaeter  that  is  capable  of  aodulating  the  intensity  of  the 
optical  beaa  when  a  voltage  or  an  electric  field  is  applied  to  an 
electro-optic  crystal .  The  structure  of  a  Mach-Zehnder 
interferoaeter  shown  in  Figure  1,  consists  of  an  input  Y-branch, 
an  output  y-branch,  and  two  parallel  waveguide  channels  that 
connect  the  two  Y-branches.  By  positioning  appropriate  electrodes 
on  or  near  the  parallel  waveguide  channels,  this  device  can  sense 
voltage  or  electric  field'.  Typically,  in  a  conventional 
electroauignetic  field  sensor,  an  external  field  to  be  aeasured  is 
picked  up  by  an  antenna  and  converted  to  a  voltage  which  is  then 
applied  to  the  electrodes  on  the  crystal.  The  electrodes  are 
connected  so  that  the  two  parallel  channels  experience  opposite 
polarities,  thereby,  leading  to  a  net  phase  shift  of  twice  the 
phase  shift  per  channel. 

In  a  variety  of  applications  that  require  the  aeasureaent  of 
electroaagnetic  fields,  the  presence  of  the  aetal  electrode  tends 
to  disturb  the  electric  field  under  aeasureaent.  In  severe  cases, 
the  close  proxiaity  of  the  electrodes  could  create  arcing,  thus 
creating  a  short  circuit.  The  electrodes  also  tend  to  Unit  the 
frequency  response  of  the  sensor  due  to  the  capacitive  nature  of 
the  electrical  circuit.  It  would,  therefore,  be  ideal  to  develop 
an  electric  field  sensor  that  does  not  require  any  aetal 
electrodes.  Such  a  device  could  be  acconplished  if  one  leg  of  the 


'C.H.  Bulner  and  W.K.  Bums,  ”  Linear  inter  feronetric 
aodulators  in  Ti:Lithiun  Niobate",  Journal  of  Lightwave  Technology 
2,  PP<  512-521  (1984). 
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Figure  1.  Ibic±i'-Zelmcler  interferoaeter  with  electrodes 


]l2UIL02>90-C-0078  8  SRIOO 

Optical  Interferoaeters  for  Sensing  BlectroMignetic  Fields 


Mach-Zehnder  were  to  be  poled  in  the  opposite  direction  to  that  of 
the  rest  of  the  crystal  and  the  waveguide.  The  device  thus  created 
would  be  aore  accurate  in  the  aeasureaent  of  the  electroaagnetic 
fields. 

We  have  developed  a  unique  method  of  creating  an  electrode-'less 
Mach-Zehnder  interferoneter  in  lithium  niobate  by  combining  a 
titanium  diffusion  step  followed  by  a  proton-exchange  process.  The 
titanium  diffusion  in  selected  sections  of  the  crystal  in  a  +Z 
oriented  substrate  is  capable  of  reverse  poling  that  section*. 
Sxibsequently,  the  Hach-Zehnder  waveguide  structure  is  created  at  a 
low  temperature  by  using  proton-exchange  waveguide  techniques* . 
This  technique  thus  avoids  the  conventional  approach  used  in  poling 
a  crystal  which  consists  of  heating  the  crystal  near  its  Curie 
temperature,  and  in  some  instances  beyond  its  Curie  temperature, 
while  applying  an  electric  field*'*.  Such  an  approach  would  be 
very  cumbersome  and  expensive. 

During  the  Phase  I  of  this  SBIR  program,  SRICO  successfully 
demonstrated  selective  domain  reversal  of  a  waveguide  channel  in  a 
+Z  oriented  LiWbO,.  Figures  2a  and  2b  show  the  and  -Y  end- 
faces  of  a  9  um  wide  Ti  diffused  channel  in  a  -t-Z  LiNbO,  substrate. 
The  end-faces  of  the  substrate  were  polished  and  etched  in  a 
solution  of  HF  and  HNO,. 

The  property  of  this  etch  solution  is  such  that  the  -fY  face  is 
etched  and  leaves  the  -Y  face  intact.  The  half-moon  shaped  region 
is  the  Ti  diffused  region.  In  Figure  2a,  the  Ti  diffused  region 
remains  smooth  whereas  the  background,  the  -fY  face,  shows  etch 
pits.  This  indicates  that  the  end-face  associated  with  the  Ti 
diffused  region  is  converted  to  -Y  orientation,  i.e.  the  Ti 
diffused  channel  is  domain  inverted  to  -Z  in  the  +Z  substrate. 


*S.  Miyazawa,  "Ferroelectric  domain  inversion  in  Ti  diffused 
lithium  niobate  optical  waveguide,"  Journal  of  Applied  Physics 
pp.  4599-4603  (1979). 

*J.L.  Jackel  et.al.,  "Proton  exchange  for  high  index 
waveguides  in  lithium  niobate,"  Topical  Meeting  on  Integrated  and 
Guided-Wave  Optics,"  Pacific  Grove,  CA,  1982. 

*K.  Nassau  et.al.,  "Ferroelectric  lithium  niobate. 
2.  Preparation  of  single  domain  crystals",  J.  Phys.  Chem.  Solids 
21,  pp.  989-996  (1966). 

*M.  Tasson  et.al.,  "Piezoelectric  study  of  poling  in  lithium 
niobate  crystals  at  temperature  close  to  the  Curie  point" , 
Ferroelectrics  12,  pp.  479-481  (1976). 
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rigiir*  2a 
<fY  End-Face 


Figure  2b 
-Y  End-Face 


Figures  2  a,b.  Ptootonicrograplis  of  +Y  and  -Y  end-faces  of  a  Ti 
diffused  +Z  oriented  LiEbO,  after  etching  in  HF:HIIO, 
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The  features  on  the  -y  end**face  are  exactly  the  opposite  as  shown 
in  Figure  2b.  In  Figure  2b,  the  Ti  diffused  region  exhibits  etch 
pits  whereas  the  background,  the  -Y  face,  remains  smooth. 


This  domain  inversion  selectively  occurs  only  in  the  Ti  diffused 
regions.  We  have  also  demonstrated  that  100%  domain  inversion  can 
be  achieved  by  choosing  appropriate  diffusion  conditions.  There 
is,  therefore,  no  degradation  in  the  electro-optic  properties  of 
the  crystal.  We  expect  the  electrode-less  device  to  have  identical 
performance  characteristics  as  the  push-pull  electroded  Nach- 
Zehnder  interferometer. 

The  specifications  for  the  electric  field  sensor  provided  in  the 
following  page  are  based  on  our  study  during  the  Phase  I  of  the 
SBIR  program.  The  performance  of  the  sensor  for  a  device  operating 
at  1,  3  and  30  6H2  are  presented  in  Table  I. 


Table  JL 


Halfwave  Voltage,  V, 

Volts 

40 

50 

67 

Substrate  Thickness 

mm 

.5 

.5 

.5 

Number  of  Segments 

1 

1 

1 

1 

10 

Segment  Interaction 
Length 

mu 

25 

20 

2.3 

Device  Bandwidth 

GHz 

DC  to  1 

DC  to  3 

30  ±  1.5 

P. 

mW 

1 

2 

4 

mW 

>2 

>4 

>8 

p 

*■  movrem 

mW 

5 

10 

10 

Wavelength 

^m 

1.3 

1.3 

1.3 

Signal  Bandwidth 

GHz 

1 

1 

1 

Dynamic  Range 

dB 

48.6 

50.83 

52.82 

E^n 

V/m 

570 

490 

460 

1.5  X  10* 

1.7  X  10* 

2.0  3?  10" 
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The  substrate  naterial  is  LiNbO,. 

The  LiNbO,  crystal  has  a  threefold  axis  of  synsetry  which  places  it 
in  the  trigonal  (rhonbohedral)  point  group  3b,  space  group  R3c. 
SoBe  of  the  physical  properties  of  LiNbO,  are  sho%m  in  Table  IT*. 


Table  II 

Properties  of  LithiuB  Miobate 


niysical  Paraaeter 

Value 

Units 

Density 

4.64 

g/CB* 

Melting  Point 

1253 

C 

Curie  Tesperature 
StoichioBetric 
Congruent 

1210 

1150 

c 

C 

Moh  Hardness 

5 

~ 

TherBal 

Conductivity 

10'“ 

cal/cB.s.c 

Pyroelectric 

Coefficient 

4  §  300K 

nC/cB*/C 

Resistivity 

10“  e  24C 

5  X  10*  §  400C 

140  f  1200C 

Ohn-CB 

Relative 

PerBittivity 

€xi/  €0  *  85  •  2 

e„/€o  -  28.7 

Cxi/eo  =  44.3 

€„/€o  =  27.9 

ClaBped 

Cleunped 

UnclaBped 

Unclasped 

*R.S.  Weis  and  T.K.  Gaylord,  "LithiuB  niobate:  SuBBary  of 
physical  properties  and  crystal  structure”.  Applied  Physics 
pp.  191-203  (1985). 
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RBFRACTIVE  I] 

The  refractive  index  of  the  crystal  depends  on  crystal  composition. 
LiNbO,  commercially  available  is  congruent  melting,  having  a  Li/Nb 
mole  ratio  of  0.946  (or  48.6  mol  %  of  LiO,) .  Also,  being  a 
negative  uniaxial  crystal,  the  extraordinary  index  n.  is  less  than 
the  ordinary  refractive  index  n,. 

The  refractive  index  values  reported  at  the  wavelength  of  633  nm 
and  at  T  =  24.5  C  are: 

n.  Stoichiometry  2.2005  n.  Congruent  Melting  2.2028 

n„  Stoichiometry  2.291  n^  Congruent  Melting  2.2866 

The  index  of  refraction  dependence  on  temperature  is  given  in  the 
paper  by  Boyd  et.al.’ 

Lithium  niobate  is  a  very  important  integrated  optic  waveguide 
material.  This  is  a  ferroelectric  material.  It  is  characterized 
by  high  piezoelectric,  electro^optlc,  pyroelectric  and 
photo-elastic  coefficients.  It  has  an  excellent  acousto-optic 
figure  of  merit.  In  this  program,  we  make  use  of  the  high 
transparency  of  this  crystal  material  in  the  wavelength  of  interest 
(visible  to  near  IR) ,  the  ability  to  create  excellent  low  loss 
waveguides  and  a  high  electro-optic  coefficient  of  the  crystal 
material . 

Below  its  Curie  temperature  of  1210  C,  the  crystal  is  in  its 
ferroelectric  phase;  and  above  T,,  it  is  in  its  paraelectric  phase. 
Since  we  will  not  be  operating  the  crystal  during  processing  or 
device  operation  beyond  its  Curie  temperature,  we  shall  limit  our 
discussion  to  its  ferroelectric  phase. 

In  its  ferroelectric  phase,  it  consists  of  planar  sheets  of  oxygen 
atoiBS  in  a  distorted  hexagonal  close-packed  configuration.  The 
octahedral  interstices  formed  in  this  structure  are  one  third 
filled  by  lithium  atoms,  one  third  filled  by  niobium  atoms,  and 
one-third  vacant.  In  the  -i-C  direction,  the  atoms  occur  in  the 

interstices  in  the  following  order:  _ ,  Mb,  vacancy,  Li,  Nb, 

vacancy,  Li,  _ .  (See  Fig.  1,  Weis  and  Gaylord). 

This  is  a  member  of  the  trigonal  crystal  system  as  it  possesses  a 
threefold  rotation  symmetry  about  its  C  axis.  It  also  exhibits  a 
mirrc.  symmetry  about  three  planes  that  are  60  degrees  apart  and 
intersect  forming  a  threefold  rotation  axis.  (See  Fig.  3,  Weis  and 
Gaylord) . 


^G.D.  Boyd  et.  al.,  "Refractive  index  as  a  function  of 
temperature  in  LiNbO,",  J.  Appl.  Phys.  pp.  1941-1943  (1967). 
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In  this  section  we  provide  a  brief  overview  of  the  physical 
principles  of  the  integrated  optic  conponents  that  Bake  up  the 
electric  field  sensor.  This  includes  a  discussion  about: 


The  integrated  optical  waveguide  that  will  be  used  to  confine 
the  optical  beam  within  the  substrate  BMterial; 

The  integrated  optical  waveguide  Y-branch  and  directional 
coupler  that  will  be  used  in  splitting  and  coBbining  a  guided 
optical  wave  propagating  in  the  waveguide; 

The  electro-optic  effect  that  will  be  used  to  change  the  index 
of  refraction  of  the  substrate  Baterial  and  thereby  affect  the 
phase  of  the  optical  beaB  propagating  in  the  optical 
waveguide; 

The  Mach-Zehnder  inter feroBeter  that  will  convert  the  phase 
shift  produced  by  the  electro-optic  effect  to  an  intensity 
variation  in  the  light  output,  thus  ened>ling  the  detection  of 
an  external  electric  field. 


OPTTCAl.  1 


JXi 


Optical  waveguides  fom  the  basic  building  block  of  integrated 
optical  devices.  The  waveguide  is  used  to  confine  the  light  beaa 
in  a  few  Bicron  layers  within  the  substrate.  A  layer  of  Baterial 
can  act  as  an  optical  waveguide  if  it  satisfies  the  following 
requireaents : 


•  The  waveguide  layer  (the  core)  and  the  substrate  Baterial  (the 
cladding)  Bust  be  transparent  at  the  wavelength  of  the  optical 
beaa; 

•  The  waveguide  layer  Bust  have  an  index  of  refraction  higher 
them  its  surrounding  Baterials; 

•  The  waveguide  layer  Bust  have  a  Biniaua  thickness  to  propagate 
at  least  one  guided  Bode. 

An  integrated  optical  waveguide  device  incorporates  a  plemar  layer 
or  a  channel  formed  in  a  material  adapted  to  a  specific 
application.  The  class  of  materials  of  interest  for  optical 
waveguide  devices  includes  glass,  silica,  or  silicon,  polymers, 
III-V  seBiconductors ,  dielectrics,  such  as  potassium  titanyl 
phosphate  (KTP)  amd  f erroelectrics ,  such  as  lithium  niobate  and 
lithium  tantalate.  Optical  waveguide  device  materials  such  as  KTP, 
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lithiua  niobate  and  lithiun  tantalate  utilize  the  electro*^ptic 
properties  of  the  crystal  saterial  for  Modulation  and  switching  of 
the  guided  beas. 

SoM  of  the  techniques  used  to  fabricate  optical  %raveguides  are 
diffusion,  depositon  and  ion-exchange. 

In  the  diffusion  process,  a  thin  layer  of  Material  is  deposited  on 
top  of  the  substrate  Material  and  diffused  into  the  Material  at 
elevated  teMperatures.  For  exanple,  an  optical  waveguide  is  forsed 
in  a  lithiuM  niobate  substrate  by  diffusing  a  layer  of  titaniun  at 
about  1025  C. 

In  the  depostion  process,  a  layer  of  high  index  Material  is 
deposited  on  top  of  the  siibstrate  Material.  The  deposition 
processes  typically  used  are  sputtering,  evaporation  and  plasMa 
deposition.  For  exasple,  a  layer  of  silicon  nitride  when  deposited 
on  top  of  an  oxidized  silicon  substrate  could  act  as  an  optical 
waveguide. 

Ion-exchange  waveguides  are  forsed  by  iMMersing  a  substrate 
Material  in  suitable  solutions  such  that  ions  in  the  solution  enter 
the  substrate  to  increase  the  refractive  index  of  the  ion-exchanged 
region.  For  exasple,  by  iamersing  soda  liMe  glass  in  a  silver 
nitrate  solution,  a  waveguide  layer  is  forsed  in  the  glass 
substrate  when  the  silver  ions  enter  the  substrate  and  replace  the 
sodiuM  ions  in  the  glass. 

There  are  two  basic  types  of  optical  waveguides.  The  first  type  is 
the  planar  waveguide  which  covers  the  entire  top  surface  of  the 
substrate  Material.  The  waveguide  is  bounded  only  in  one  direction 
by  the  thickness  of  the  waveguide  layer.  The  optical  beaM  is 
capable  of  propagating  in  any  direction  in  the  plane  of  the 
waveguide.  This  planar  waveguide  geoMetry  is  used  in  acousto-optic 
devices  in  idiich  the  guided  wave  is  diffracted  at  Many  different 
angles  by  the  acoustic  surface  waves*. 

The  second  type  of  optical  waveguide  is  the  channel  waveguide  in 
%rhich  the  light  beaM  is  confined  to  propagate  in  a  specific  path 
defined  by  the  channel.  The  waveguide  is  bounded  in  both  the 
thickness  and  the  width,  and  the  light  beaM  propagates  along  the 
length  of  the  channel,  ^e  light  beu  is  thus  well  confined  within 
the  channel  region.  The  channels  are  fabricated  by  confining  the 
depostion  or  diffusion  process  to  the  required  width  by  using  the 
photolithograpic  process.  Cheumel  waveguides  are  used  to  direct 
the  optical  beaM  to  specific  locations  by  forsing  the  desired 


*C.S.  Tsai,  et.  al.,  "Hide-band  guided-%mve  acousto-optic 
Bragg  diffraction  emd  devices  using  Multiple  tilted  surface 
acoustic  waves",  Proc.  IEEE  318  (1976). 
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optical  circuits  on  the  substrate.  The  electric  field  sensor  under 
developaent  utilises  channel  waveguides  in  a  lithius  nic^te 
electro-optic  crystal  saterial. 


The  Hach-Zehnder  interferoseter  uses  a  Y-branch  at  its  input  and 
output.  At  the  input,  the  Y-kuranch  acts  as  a  power  divider  to 
excite  the  two  legs  of  the  interfercsMter.  For  a  perfectly 
sysMtric  junction,  the  power  will  be  divided  equally  between  the 
two  eras  as  shown  in  Figure  3a.  Any  departure  froa  syaaetry  will 
result  in  unequal  power  splitting. 

At  the  output  of  the  Hach-Zehnder,  the  Y-branch  acts  to  coabine  the 
light  propagating  in  the  t%fo  legs  of  the  interferoaeter.  For  equal 
inputs  into  the  tiro  aras,  the  output  will  vary  froa  sero  to  the  sua 
of  the  two  input  intensities,  depending  on  the  relative  phase  of 
the  two  inputs.  If  the  two  inputs  are  equal  and  in  phase,  the 
output  of  the  y-branch  will  be  equal  to  the  sua  of  the  two  input 
intensities  as  shown  in  Figure  3b.  But,  if  the  am>litude  of  the 
two  inputs  are  equal  but  opposite  in  phase,  the  output  of  the  Y- 
brcuich  will  be  zero.  The  energy  that  does  not  reach  the  output  ara 
is  radiated  into  the  substrate  as  shoim  in  Figure  3c.  In 
principle,  the  Y-branch  is  a  four  port  device  —  the  three  aras  of 
the  y-branch  account  for  three  ports  and  the  substrate  is  the 
fourth  port. 


Another  waveguide  coaponent  that  could  be  used  as  a  power  splitter 
or  a  coabiner  is  a  directional  coupler  shown  in  Figure  4 .  A 
directional  coupler  consists  of  two  parallel  channel  waveguides 
fabricated  close  together  such  that  the  two  waveguides  interact 
with  each  other.  The  evanescent  optical  fields  in  the  two 
waveguides  overlap  such  that  the  energy  is  transferred  between  each 
other.  The  coupling  of  light  between  the  two  traveguides  depend  on: 

•  The  separation  between  the  two  waveguides; 

•  The  length  over  which  interaction  takes  place  between  the  two 
channel  waveguides; 

•  The  index  of  refraction  of  the  waveguide  and  the  surrounding 
■ediuB. 

The  directional  coupler  is  a  four  port  device.  One  sajor 
difference  between  the  directional  coupler  and  the  Y-branch  is  that 
all  four  ports  of  the  directional  coupler  are  waveguides.  Hence, 
there  is  no  energy  wasted  in  the  susbtrate  due  to  radiation. 
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Th«  directional  coupler  could  therefore  be  used  to  replace  the 
output  y-'branch  to  reconbine  the  optical  signal  without  any  loss  of 
energy.  The  directional  coupler  could  be  designed  and  feUt>ricated 
such  that  the  t%io  output  ports  have  equal  intensities.  One 
additicNial  feature  of  the  directional  coupler  is  that  the  ttio 
outputs  are  always  out  of  phase.  We  will  later  aake  use  of  this 
property  of  the  directional  coupler  to  increase  the  carrier-to~ 
noise  ratio  of  the  electric  field  sensor. 


The  key  paraaeter  of  interest  for  this  interferoaeter  device  is  the 
linear  electro-optic  effect,  also  known  as  the  Pockels  Effect.  The 
linear  electro-optic  effect  is  aanifested  in  a  non-centro  syaaetric 
crystal  «dien  an  electric  field  is  present.  The  relationship 
between  the  refractive  index  change  and  the  electric  field  is 

^  -  E.  (1) 


where,  A(l/n^)i^  is  the  second-rank  tensor  describing  the  change  in 
relative  pemittivity,  E^  is  the  k-th  coaponent  of  the  electric 
field  vector  and  r^^  is  the  third-rank  linear  electro-optic 
coefficient  tensor.  The  electro-optic  tensor  of  27  eleaents  can  be 
reduced  to  a  6x3  aatrix.  Lithiua  niobate  has  only  four  independent 
coefficients  (r.,,  r,,,  rx,  and  r,,).  The  ri^^  tensor  in  reduced 
subscript  notation  is  given  by: 


rax=  6.8  X  10-”  a/V,  rx,=  10  x  10'”  a/v, 

r,x=  30.8  X  10"”  a/v,  r^*  32.6  x  10-“  a/V. 

These  values  are  at  the  wavelength  of  633  na. 
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A  symetric  Mach-Zehnder  (NZ)  inter feroaeter  device  is  shown  in 
Figure  1.  This  is  an  electroded  device.  This  device,  when  driven 
by  a  dipole  antenna,  would  function  as  an  electric  field  sensor. 
The  dipole  antenna  converts  the  electric  field  to  a  voltage  output. 
This  voltage  is  applied  to  the  electrodes  on  the  device.  The 
electrode  arrangement  shown  in  Figure  1  is  for  a  -Z  oriented  LiNbO, 
crystal . 

The  waveguide  structure  shown  consists  of  an  input  channel  that 
propagates  the  light  into  the  input  Y-branch.  At  the  Y-branch,  the 
light  is  equally  split  between  the  two  channels.  The  two  channels 
branch  out  at  an  angle  of  one  degree  to  a  separation  of  50  ^m.  The 
two  channels  are  bent  so  that  they  are  parallel  to  each  other  over 
the  length  L.  The  two  parallel  channels  then  recombine  at  the 
output  Y~branch  and  transmit  the  recombined  light  into  the  output 
port.  In  the  symmetric  device  shown,  the  optical  path  lengths 
associated  with  the  light  in  the  upper  and  lower  half  of  the  MZ 
interferometer  are  equal.  In  the  absence  of  an  external  voltage, 
the  two  optical  beams  experience  equal  phase  shift,  leading  to  a 
net  phase  shift  of  zero. 

The  electrodes  are  placed  such  that  the  electric  field  generated  in 
the  crystal  due  to  the  applied  voltage  is  in  the  Z  direction.  The 
electric  field  is  oppositely  directed  in  the  two  legs  of  the 
interferometer . 

The  effect  of  this  electric  field  is  to  change  the  refractive  index 
of  the  material  through  the  linear  electro-optic  effect.  Since  the 
electric  field  is  only  in  the  Z  direction,  E.,  Ey  -  0.  Only  E,  is 
a  non-zero  component  of  the  electric  field  vector.  The  index 
change  along  each  direction  can  be  calculated  using  the  electro¬ 
optic  tensor  r^^.  If  we  let  1,2,3  correspond .  to  x,y,z  co¬ 
ordinates,  we  have 
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Based  on  the  above  equations,  the  refractive  index  of  the  crystal 
due  to  the  application  of  an  external  electric  field  can  be 
calculated. 


In  the  case  of  LiNbO,,  the  refractive  index  along  the  x  and  y 
directions  are  the  seuie,  referred  to  as  the  ordinary  index,  n*. 
The  refractive  index  along  the  z  axis  (or  the  c  axis)  ,  referred  to 
as  the  extraordinary  index,  is  n,.  If  we  restrict  the  light 
propagating  in  the  crystal  to  the  TN  sode,  which  has  its  E  vector 
in  the  z  direction,  we  can  then  neglect  the  effect  of  increase  in 
the  ordinary  index,  n,.  The  extraordinary  index  due  to  an  electric 
field,  E,,  is  given  by 


+  ^33 


(5) 


,  1  *  ^33 

"•Lo 


(6) 


^33 


(7) 


1 

2 


"i  ^33 


(«) 


This  is  the  change  in  the  refractive  index  in  each  leg  of  the 
interferometer  due  to  an  external  electric  field. 

If  V  is  the  voltage  applied  to  the  electrodes  and  d  is  the  gap 
between  them,  the  electric  field  E.  is  given  by 


Then, 


An. 


1 

2 


n: 


(10) 


The  phase  change  in  the  optical  beam  propagating  in  each  leg  of  the 
interferometer  is 
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=  --J  V  ^33  -g  ^ 


(11) 


where,  X  is  the  wavelength  of  the  light  propagating  in  the 
waveguide. 

The  device  configuration  8ho%m  in  Figure  1  operates  in  a  push-pull 
fashion,  ^e  phase  chemge  associate  with  the  second  leg  of  the 
interferoeeter  is  opposite  in  sign  to  that  of  the  first  leg,  i.e., 

A*,  =  (12) 

The  total  phase  change  is 

|A*.  -  A*a|  -  |A*,| 

-  2  ^  n2  r„  i 


The  optical  signal  output  of  the  interferometer  is  given  by 

p.  .  cos>(-!^)  (I*) 


(13) 

(14) 


A  plot  of  this  equation  is  shown  in  Figure  5. 

The  device  operating  point  is  at  A#«  »  0.  This  is  a  non-linear 
portion  of  the  cosine  squared  curve.  The  output  signal  is  not  very 
sensitive  to  phase  change  at  this  operating  point.  Higher 
sensitivity  is  obtained  by  operating  the  device  in  the  linear 
regime,  ffe  observe  from  Figure  5  that  the  linear  portion  of  the 
curve  occurs  for  a  phase  change  of  w/2  radian  where  the  output 
signal  is  one  half  its  maximum. 

By  operating  the  device  at  this  point,  the  device  operates  in  a 
linear  fashion.  Since  the  basic  curve  is  a  cosine  squared 
function,  linearity  is  maintained  only  over  a  small  {diase 
deviation.  The  maximum  phase  deviation  depends  on  the  harmonic 
distortion  allowed  in  the  device  performance. 
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DOMAIM  IMVERSIOII 


Ferroelectric  lithiim  niobate  exhibits  spontaneous  polarization 
below  its  Curie  temperature,  T^.  The  orientation  of  this 
spontaneous  polarization,  P.,  can  be  affected,  and  in  many 
situations  reversed,  by  the  application  of  an  external  electric 
field,  a  temperature  gradient,  a  concentration  gradient,  or  an 
electron  beam. 

While  all  four  techniques  have  shown  success  in  domain  reversal  in 
lithium  niobate,  it  is  only  the  concentration  gradient  technique 
that  lends  itself  easily  for  applications  in  integrated  optics.  In 
integrated  optics,  selective  sections  of  an  optical  circuit  would 
require  that  the  orientation  of  the  spontaneous  polarization  be 
reversed  compared  to  the  bulk  crystal  material.  This  is  easily 
accomplished  by  the  concentration  gradient  effect  by  the  mechanism 
of  high  temperature  diffusion  of  dopants  into  the  crystal.  While 
elaborate  schemes  can  be  designed  to  apply  an  electric  field  or 
temperature  gradient  or  an  electron  beeun  in  selective  regions  of 
the  crystal,  this  is  quite  cumbersome  and  does  not  lead  to  a 
manufactured)! e  process. 

We  will  provide  a  brief  e}cplanation  of  the  electric  field,  the 
temperature  gradient  and  the  electron  beam  effects  and  explain  in 
detail  the  effects  of  the  concentration  gradient  due  to  doping. 


£LB£ZBI£_E1£UI 


A  crystal  of  LiNbO,  is  heated  above  its  Curie  point.  During  the 
cooling  process,  an  electric  field  of  1  V/cm  is  applied.  When  a 
temperature  T.  is  reached,  such  that  T^  is  slightly  greater  than  T^, 
the  electric  field  is  turned  off  and  the  crystal  continues  to  cool 
to  room  temperature.  In  this  situation,  the  crystal  has  a 
spontaneous  polarization  that  is  opposite  to  that  of  the  starting 
material.  Domain  inversion  is  thus  achieved  in  the  crystal. 


TEMPERATDRE  GRADIENT  EFFECT 

A  temperature  gradient  of  100  C/cm  across  the  sample  thickness 
results  in  a  similar  effect  to  those  observed  under  electric 
fields.  When  a  sample  is  heated  above  its  Curie  temperature  and 
cooled  while  maintaining  a  temperature  gradient,  the  material  gets 
poled  in  a  sense  which  is  related  to  the  thermal  gradient  sign.  If 
the  same  is  now  quickly  heated  above  T^  and  cooled  to  room 
temperature  without  any  temperature  gradient,  the  sample  gets  poled 
in  a  sense  which  is  opposite  to  that  of  the  starting  material.  We 
have  thus  achieved  domain  inversion  in  the  crystal. 
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DQMAIM  REVERSAL  OSlllG  ELECTROW  BEAM 


This  technique  was  published  by  Hayccx:k  and  Townsend*  in  1986.  It 
permits  domain  reversal  on  both  the  +C  and  -C  surface  of  LiNbO,  and 
LiTaO,  at  temperatures  much  below  their  Curie  temperatures.  This 
also  requires  simultaneous  application  of  an  electric  field  during 
the  electron  bombardment  of  the  crystal. 

The  domain  reversal  was  achieved  in  LiNbO,  at  600  C  (To  =  1140  C) 
with  an  applied  electric  field  of  10  V/cm  and  a  beam  energy  of 
1.8  MeV  for  a  1  mm  thick  crystal. 

The  domain  reversal  was  achieved  in  LiTaO,  at  400  C  (T^  «  620  C) 
with  an  applied  electric  field  of  900  V/cm  and  a  beam  energy  of 
1.8  NeV  for  a  1  mm  thick  crystal. 

This  approach  would  allow  domain  reversal  in  selected  regions  of  a 
crystal  to  form  optical  waveguide  circuits.  This  is  possible  since 
the  crystal  is  covered  by  a  gold  film  patterned  by  photolithography 
to  protect  the  rest  of  the  crystal  from  being  bombarded  by  the 
electron  beam.  An  electric  field  is  applied  to  this  gold  electiode 
during  the  process. 

VRiile  this  technique  offers  the  choice  of  -<-C  and  -C  surfaces  and  a 
lower  temperature  processing,  it  requires  access  to  an  electron 
beam  source.  Using  the  e-beam  process.  Keys  et.  al.***  have  created 
a  grating  structure  for  obtaining  SHG  in  LiNbO,. 


COMCKHTRATION  GRADIENT 


Concentration  gradient  is  created  in  bulk  crystals  by  introducing 
dopants  in  the  crystal  growth  process.  For  integrated  optic 
devices,  dopants  are  introduced  into  the  crystal  in  selected 
waveguide  regions  by  diffusion.  Niyazawa  discovered  domain 
inversion  in  LiNbO,  when  Ti  metal  was  diffused  into  the  crystal. 
In  particular,  it  is  only  the  diffusion  of  Ti  into  the  +C  surface 
layer  that  resulted  in  domain  inversion.  The  -C  surface  showed 
almost  no  domain  inversion  due  to  Ti  diffusion.  He  was  able  to 
create  both  slab  waveguide  and  stripe  waveguide  in  the  +c  surface 
of  lithium  niobate.  In  a  slab  waveguide  configuration,  a  layer  of 
Ti  is  diffused  on  the  complete  surface  of  the  crystal.  In  this 
case,  light  confinement  is  possible  only  in  the  waveguide 


*P.W.  Haycock  and  P.D.  Townsend,  "A  method  of  poling  LiNbO,  and 
LiTaO,  below  T*,",  Applied  Phys.  Lett.  A&,  pp.  698-700  (1986). 

'*11. W.  Keys  et.  al.,  "Fabrication  of  domain  reversed  gratings 
for  SHG  in  LiNbO,  by  electron  beam  bombardment".  Electronics  Lett. 
2£/PP*  188-189  (1990). 
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thickness.  In  a  stripe  waveguide  configuration,  Ti  is  patterned  to 
a  width  of  a  few  microns  and  then  diffused.  In  this  case,  light 
confinement  is  possible  in  both  the  lateral  and  vertical 
directions. 

A  simple  chemical  etch  technique  was  used  to  reveal  the  domain 
structure  in  the  Ti -diffused  layers.  A  boiling  solution,  110  C,  of 
HNO,  and  HF  (2:1)  was  used  to  etch  the  sample  over  three  minutes. 

An  undiffused  crystal  of  LiNbO,  as  received  from  the  manufacturer 
will  exhibit  deeply  etched  patterns  on  the  -C  surface  and  scarcely 
etched  pattern  on  the  -fC  surface.  After  Ti  diffusion,  the  -t-C 
surface  shows  deeply  etched  pattern.  The  etch  pattern  on  the  Ti 
diffused  -C  surface  is  no  different  from  the  original  crystal 
surface.  This  indicates  that  the  Ti  diffused  layer  behaves  in  the 
fashion  of  a  -C  surface  leading  to  the  conclusion  that  a  domain 
inversion  has  taken  place  in  that  region. 

The  crystal  was  wrapped  in  a  Pt  foil  to  suppress  out  diffusion  of 
Li  at  high  temperatures.  A  slab  waveguide  created  by  diffusing 
500A  Ti  at  1050  C  for  10  hours  could  support  three  (3)  modes  at  the 
wavelength  of  the  He-Ne  laser  (A.  >>  6328A). 

An  experiment  to  characterize  the  degree  of  depoling  (or  domain 
inversion)  was  carried  out  by  diffusing  Ti  at  several  temperatures 
from  950  C  to  1100  C  for  5  h  and  10  h.  Only  partial  depoling 
occurred  below  1020  C.  For  the  shorter  diffusion  time  of  5  h,  100% 
depoling  occurred  at  lower  temperature  compared  to  10  h  diffusion 
time.  At  temperatures  beyond  1020  C,  100%  depoling  occurred  for 
both  5  and  10  h  diffusion  time.  The  average  diffusion  depth  was 
found  to  be  10  microns.  The  Ti  concentration  at  the  top  surface  is 
higher  for  a  5  h  diffusion  than  for  a  10  h  diffusion  process  at  the 
same  temperature.  For  temperatures  where  100  %  depoling  does  not 
take  place,  it  was  observed  that  the  shorter  diffusion  time  process 
yields  a  higher  degree  of  depoling,  leading  to  the  conclusion  that 
the  domain  inversion  is  related  to  the  Ti  concentration. 

Stripe  waveguides  were  also  created  for  several  widths  ranging  from 
4  to  10  ^m.  It  was  noticed  that  4  and  5  ixm  wide  guides  exhibited 
partial  domain  inversion  whereas  guide  widths  of  6  -  10  |im  were 
100%  depoled.  O^e  Ti  concentration  in  the  4  and  5  |im  wide  guides 
are  not  believed  to  be  sufficient  enough  to  create  100%  depoling. 
Hence,  a  thicker  Ti  layer  would  achieve  the  desired  100%  domain 
inversion  in  these  4  and  5  (im  guide  widths. 
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The  process  requires  defining  Ti  aetal  in  selected  regions  of  the 
substrate  for  subsequent  diffusion  operation.  Standard 
photolithographic  process  techniques  developed  in  aicroelectronic 
device  technology  is  used  for  Ti  pattern  definition. 

In  the  process  shown  in  Figure  6,  a  layer  of  photoresist  is  spin 
coated  on  to  the  substrate.  A  photomask  that  has  the  required 
pattern  «rritten  on  it  is  placed  on  top  of  this  substrate  in  a  mask 
aligner.  The  photoresist  is  exposed  to  a  UV  light  source  through 
the  photomask.  When  the  substrate  is  processed  in  developer 
solution,  the  photoresist  is  etched  off  from  regions  that  were 
exposed  through  the  photomask.  A  layer  of  Ti,  500  A  to  1000  A,  is 
deposited  by  evaporation  to  cover  the  entire  siurface  of  the 
substrate.  When  this  substrate  is  soaked  in  acetone,  the 
photoresist  is  removed  and  thus  lifts  off  the  Ti  metal  that  is  on 
top  of  the  photoresist.  Ti  metal  is  present  only  in  regions  of  the 
substrate  that  had  the  photoresist  etched  in  the  earlier  step.  The 
substrate  is  cleaned  and  is  ready  for  diffusion. 

The  diffusion  process  is  carried  out  in  an  open  tube  furnace  at 
around  1025  C.  The  substrate  is  placed  on  an  aliunina  boat  which  is 
positioned  in  the  center  zone  of  the  ftimace.  The  furnace 
atmosphere  is  oxygen  saturated  with  water  vapor  to  prevent 
outdif fusion  of  titanium.  Typical  oxygen  flow  rate  is  around  1  to 
1.5  SLM  (standard  liters  per  minute)  and  the  water  temperature  is 
greater  than  60  C. 

The  furnace  temperature  is  increased  gradually  to  the  desired 
diffusion  temperature  of  1025  C  over  a  2  hour  period  so  as  not  to 
thermally  shock  the  crystal.  At  the  end  of  the  diffusion  period, 
the  furnace  is  turned  off  and  allowed  to  cool  to  room  temperature 
before  retrieving  the  substrate  from  the  ftimace. 
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Spin  Photoresist 


Lithium  Niobate 


Photolithography 

Expose  and  Develop  Photoresist 


Deposit  Titanium 

500  to  1000  Angstroms 


Lift-Off  Process 
SoakinAoetorw 


Diffuse  Titanium  for  Domain  Reversal 
1025  C,  8  Hours 


Figure  6.  Process  sequence  for  creating  dosain  inversicm 
in  a  -fZ  oriented  lithius  nicfoate  substrate 
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Feasibility  experisents  were  carried  out  during  the  Phase  I  SBIR 
prograe  to  deeonstrate  selective  polarization  reversal  in  a  -fZ 
oriented  LiNbO,  substrate.  The  process  conditions  were  deterained 
to  obtain  100%  dozuiin  reversal  in  the  Ti  diffused  regions. 


The  experiments  to  characterize  the  degree  of  depoling  (or  doaain 
inversion)  were  carried  out  by  diffusing  Ti  in  LiNbO,  at  a 
teaperature  beyond  1000  C.  A  simple  chemical  etch  technique  was 
used  to  reveal  the  domain  structure  in  the  Ti -diffused  layers.  A 
solution  of  HNO,  «md  HF  (2:1)  at  a  temperature  of  80  to  90  C  was 
used  to  etch  the  sample  over  three  minutes. 

First,  a  -Z  oriented  LiNbO,  substrate  was  diffused  with  lOOOA  of  Ti 
at  1025  C  for  8  hours.  After  diffusion,  the  crystal  was  etched  in 
the  above  mentioned  solution.  The  etch  pits  thus  produced,  shown 
in  Figure  7,  are  typical  of  the  -Z  oriented  LiNbO,  crystal 
surface.  We  will  use  Figure  7  as  a  reference  to  measure  the  degree 
of  domain  inversion  in  the  +Z  oriented  crystal. 


Figure  7.  ~Z  oriented  lithium  ni<^te  crystal  diffused  trith  lOOOA 
of  titanium  and  subsequently  etched  in  HF:HMO,  at  80  to 
90  C  for  3  minutes 
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The  +Z  oriented  LiNbO,  is  normally  resistant  to  the  acid  etch 
solution  described  above.  But,  after  Ti  diffusion,  the  -fZ  surface 
exhibits  deeply  etched  patterns.  Figures  8  and  9  show  the  etch 
patterns  of  a  Ti  diffused  +Z  LiNbO,  surface.  The  thickness  of  the 
Ti  were  1000  and  SOOA  respectively  and  the  diffusion  was  carried 
out  at  1025  C  for  8  hours.  The  etch  pattern  on  the  Ti  diffused  +Z 
surface.  Figure  8,  is  identical  to  that  of  the  -Z  surface  as  shown 
in  Figure  7.  This  indicates  that  the  Ti  diffused  layer  in  the  -fZ 
surface  behaves  in  the  fashion  of  a  -Z  surface  leading  to  the 
conclusion  that  a  domain  inversion  has  taken  place  in  that  region. 

The  degree  of  domain  inversion  is  estimated  from  the  etch  pit 
density  shown  in  Figures  8  and  9.  Since  the  etch  pits  of  Ficmres  7 
and  8  are  identical,  we  conclude  that  the  diffusion  of  lOOOA  of  Ti 
into  a  +Z  surface  has  led  to  a  100%  domain  reversal  to  a  -Z 
oriented  surface.  The  etch  pits  shown  in  Figure  9  associated  with 
500A  of  Ti  diffused  into  a  -fZ  surface  exhibit  a  much  lower  etch  pit 
density  compared  to  Figures  7  and  8,  indicating  only  a  partial 
domain  inversion  for  this  procer^  condition. 


Figure  8.  +Z  LiNbO, 
lOOOA  Ti  0  1025  C,  8h 
Etc>*ed  in  HF:HNO, 


Figure  9.  +Z  LiNbO, 
500A  Ti  0  1025  C,  8h 
Etched  in  HFxEDK), 
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The  diffusion  operation  was  also  carried  out  at  a  teaperature  of 
1050  C  for  8  hours  for  Ti  thickness  of  lOOOA.  The  higher  diffusion 
tenperature  leads  to  further  diffusion  of  Ti  and,  thus,  a  decrease 
in  the  concentration  gradient.  This  reduction  in  the  titaniua 
concentration  gradient  leads  to  a  less  than  100%  dosain  inversion. 
The  etch  pits  produced  in  this  sasple  after  etching  in  HF:HNO,  is 
sho%m  in  Figure  10.  The  etch  pit  density  is  significantly  lower 
than  that  of  the  reference  sample  shown  in  Figure  7.  The  degree  of 
depoling  in  this  sample  is  estimated  to  be  less  than  50%. 


Figure  10.  +z  oriented  lithium  niobate  crystal  diffused 
with  lOOOA  titaniua  at  1050  C  for  8  hours  and 
subsequently  etched  in  HP;BIiO,  at  80  to  90  C 
for  3  minutes 


Titanium  stripes  were  also  created  for  several  widths  remging  from 
2  to  15  ^ra.  The  thickness  of  the  Ti  was  lOOOA,  and  the  diffusion 
was  carried  out  at  1025  C  for  8  hours.  It  was  noticed  that  widths 
smaller  than  5  microns  exhibited  peurtial  domain  inversion  whereas 
widths  6  microns  and  greater  were  100%  depoled.  Figure  11  shows 
the  etch  pits  produced  in  a  9  ^m  wide  Ti  diffused  channel  in  a  +Z 
oriented  surface.  The  density  of  etch  pits  are  similar  to  that 
shown  in  Figures  7  and  8,  indicating  100%  domain  inversion  in  the 
Ti  diffused  stripe  region.  Also,  note  that  the  surface  outside  the 
Ti  diffused  region  does  not  exhibit  etch  pits.  The  non-Ti  diffused 
region  merely  exhibits  scratches  which  is  a  result  of  the  acid  etch 
process .  The  cross  sections  of  the  same  seunple  were  shown  in 
Figures  2a  and  2b.  The  9  pun  wide  stripe  increases  to  a  width  of 
18  ^m  due  to  lateral  diffusion.  Since  the  two  legs  of  the  Mach 
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Zehiid«r  interferometer  are  separated  by  at  least  50  the  lateral 
diffusion  would  have  no  effect  on  the  device  performance. 

we  conclude  that  selective  100%  domain  reversal  is  accomplished  in 
a  +Z  LiNbO,  surface  by  diffusing  a  lOOOA  thick  Ti  at  1025  C  for 
8  hours.  The  width  of  the  Ti  stripe  should  at  least  be  6  microns. 


Figure  11.  +Z  oriented  lithium  niobate  crystal  showing 
9  micron  wide  channel  diffused  with  lOOOA  titanium 
at  1025  C  for  8  hours  and  subsequently  etched  in 
HF:HIiO,  at  80  to  90  C  for  3  minutes 
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-KXCHAIIGB  MAVEGOIDB  PROCESS 

The  Mach~Zehnder  interferoeeter  waveguide  structure  is  to  be 
fabricated  in  a  +2  oriented  LiNbO,  substrate.  Selective  sections 
of  this  crystal  are  reverse  poled  using  the  high  teeperature 
diffusion  process  explained  earlier.  One  leg  of  the  interferoaeter 
is  aligned  to  be  within  this  reverse  poled  section  of  the  crystal. 

The  process  of  creating  the  waveguide  structure  Bust  be  such  that 
it  does  not  degrade  the  previous  process  performed  on  the  crystal. 
There  are  two  well  established  processes  for  feUbricating  waveguides 
in  LiNbO,.  These  are  high  temperature  diffusion  of  titanium  and 
proton  exchange  in  organic  acids. 

In  the  specific  crystal  orientation  chosen  for  this  device,  the 
process  of  Ti  diffusion  leads  to  domain  inversion  for  temperatures 
higher  than  1025  C.  It  is  reported  that  at  1000  C  diffusion 
temperature,  no  domain  inversion  takes  place.  We  could  therefore 
fabricate  waveguides  at  this  teiqperature.  Unfortunately,  the 
diffusion  time  required  would  be  very  high,  greater  than  10  hours 
to  obtain  waveguides  that  operate  at  1.3  iim  wavelength.  This 
second  diffusion  step  would  further  diffuse  the  Ti  already  present 
in  the  domain  inverted  region  and  could  severely  degrade  the 
waveguide  performance  in  this  region.  At  the  end  of  this  two<-step 
diffusion  process,  the  crystal  would  have  tundergone  a  total  of  18 
hours  of  high  temperature  processing  >  8  hours  at  1025  C  and  10 
hours  at  1000  C.  The  overall  substrate  could  be  affected  by  such 
a  long  exposure  to  high  temperature.  For  example,  this  could  lead 
to  excessive  outdif fusion  of  Li  from  the  LiNbO,  leading  to  unwanted 
surface  guides. 

The  second  approach  to  waveguide  fabrication  is  the  so  called 
proton-exchange  process.  The  waveguide  process  temperature  is  in 
the  range  of  200  to  400  C  and  would  therefore  not  affect  the  Ti 
diffused  domain  inverted  regions  or  the  crystal  substrate.  Protons 
from  benzoic  acid  exchange  with  Li  in  LiNbO,  to  form  the  waveguide. 
This  was  first  described  by  J.L.  Jackel  of  Bell  Laboratories.  This 
is  a  low  temperature  process  that  allows  for  simple  waveguide 
fed>rication  techniques.  In  addition,  there  are  two  other  distinct 
advimtages  to  using  this  technique. 

First,  the  waveguide  thus  fabricated  is  truly  single  mode.  This 
process  only  increases  the  extraordinary  refractive  index,  n.,  with 
a  slight  decrease  in  the  index  of  the  ordinary  ray,  n,.  This 
permits  only  a  single  polarization  to  propagate  in  the  waveguide 
which  is  of  importance  in  a  Nach-Zehnder  interferometer,  this  being 
a  polarization  sensitive  device.  Second,  the  proton-exchange 
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waveguides  exhibit  higher  optical  dasage  threshold^^  cc»pared  to 
Ti  diffused  waveguides.  Hence,  this  would  allow  greater  optical 
power  to  propagate  in  the  proton  exchanged  waveguide,  thus 
resulting  in  greater  device  sensitivity. 

We,  therefore,  select  the  proton-exchange  process  to  fabricate 
waveguides  for  this  Mach-Zehnder  interferometer  device. 

The  proton-exchange  (PE)  process  in  LiNbO,  can  be  broadly  described 
as  the  depletion  of  lithium  ions  and  its  replacement  with  hydrogen 
ions  in  the  crystal.  This  occtirs  when  the  crystal  is  immersed  in 
a  protonic  source  such  as  benzoic  acid. 

In  order  to  obtain  optical  waveguides,  only  a  partial  exchai^e  is 
required.  In  fact,  100%  exchange  of  H*  for  Li*^  leads  to  changes 
in  the  crystal  structure  and  cracking  of  the  crystal  due  to 
excessive  stress. 

The  reaction  mechanism  is 

LiiabO,  ♦  (CfH^COOM)  ^  (C^H^COOLi)  ^  (16) 


Successful  waveguide  formation  is  achieved  with  0.8  ^  x  St  0. 

The  exchange  process  is  accomplished  by  using  either  pure  benzoic 
acid  or  benzoic  acid  diluted  with  lithium  benzoate.  The  process 
temperatiure  is  in  the  range  of  200  to  250  C.  The  process  time 
depends  on  the  type  of  wavciguide  i.e.  single  or  multimode  and  its 
wavelength  of  operation.  The  exchange  time  is  typically  30  minutes 
to  several  hours,  again  depending  on  the  type  of  waveguide 
required. 

This  PE  waveguide  is  known  to  be  unstable  due  to  further  migration 
of  protons  within  the  thickness  of  the  guiding  layer.  The 
refractive  index  profile  has  also  been  observed  to  vary  with 
tine”.  The  waveguides  fabricated  in  pure  benzoic  acid  exhibit 
increased  propagation  loss  and  scattering  loss.  Earlier  devices 
fabricated  by  the  PE  process  also  exhibited  a  degradation  in  the 
electro-optic  effect. 


”K.K.  Wong,  "Electro-optic  frequency  translators  in  lithiiu 
niobate  f2d>ricated  by  titanium  indiffusion  and  proton  exchange", 
Ph.D.  Thesis,  University  of  Glasgow,  Glasgow,  Scotleuid,  1989. 

^%.K.  Wong  et.  al.,  "Characterization  of  proton  exchange  slab 
optical  waveguides  in  X-cut  LiNbo,",  lEE  Proceddings,  Part  J  111, 
pp.  113-117  (1986). 


DhAL02-90-C-0078  34  SRIOO 

Optical  Interferometers  for  Sensing  Electromagnetic  Fields 


The  above  problees  were  overcoee  by  Suchoski  et.  al.^*  by 
perforaing  an  annealing  operation  after  the  proton-exchange 
process.  The  annealing  operation  vas  carried  out  in  flowing  oxygen 
at  350  C  for  four  hours  to  obtain  waveguides  at  1.55  im  wavelength. 
Ibis  draaatically  reduced  the  waveguide  loss,  was  stable  over  a 
teaperature  range  of  25  -  150  C,  and  did  not  exhibit  any 
degradation  in  the  electro-optic  effect. 


PBQCIga  SBQBgaCB 

The  process  sequence  for  an  electrode-less  Nach-Zehnder 
interferoaeter  is  shown  in  Figure  12.  The  LiNbO,  substrate  has 
a  doaain  inverted  Ti  diffused  region  as  shown  in  this  figure.  A 
thin  layer  of  aluainua,  1000  A,  is  deposited  on  the  surface  of  the 
crystal.  Next,  a  layer  of  photoresist  is  coated  and  the  waveguide 
pattern  is  defined  on  the  photoresist  using  i^otolithography.  In 
the  photolithographic  step,  a  photoaask,  %diich  has  the  waveguide 
pattern  written  on  it,  is  placed  on  top  of  the  photoresist  coated 
substrate.  nie  photoresist  is  exposed,  through  the  photoaask, 
using  a  UV  light  source.  After  exposure,  the  substrate  is 
developed  in  a  developer  solution.  The  j^otoresist  is  thus  reaoved 
in  areas  tdiere  the  waveguide  is  to  be  defined.  The  underlying 
aluainua  layer  is  etched  using  photoresist  as  a  protective  layer  to 
prevent  the  aluainua  froa  etching  in  areas  underneath  the 
photoresist.  The  substrate  is  then  iaaersed  in  a  benzoic  acid  aelt 
at  200  C.  The  proton  exchange  takes  place  only  in  regions  where 
the  aluainua  has  been  etched.  The  aluainua  layer  prevents  the  PE 
process  taking  place  underneath  it.  The  aluainua  layer  is 
coapletely  reaoved  in  phosphoric  acid.  The  crystal  is  then 
emnealed  in  flowing  oxygen  at  350  C  for  4  hours  to  fora  the  desired 
waveguide  device. 

The  overall  process  to  create  an  electrode-less  Nach-Zehnder 
interferoaeter  electric  field  sensor  consists  of  a  high  teaperature 
Ti  diffusion  process  followed  by  a  proton  exchange  process.  A 
sketch  of  this  device  structure  is  shown  in  Figure  13.  The 
proton-  exchange  waveguide  in  the  lower  leg  of  the  interferoaeter 
is  placed  within  the  Ti  diffused  doaain  inverted  region. 

In  the  device  shown,  there  is  an  ed>rupt  change  in  the  propagation 
constant  of  the  optical  wave  as  it  enters  and  leaves  the  Ti 
diffused  region.  The  propagation  constant  of  the  purely  PE  region 
could  be  different  than  that  of  the  PE  region  within  the  Ti 
diffused  region  in  the  range  of  5xl0~*/|in.  This  could  lead  to 
increased  light  scattering  and  loss  at  the  junction  region. 


^’P.6.  Suchoski  Jr. ,  "Passive  polarization/polarization 
maintaining  components  based  on  proton-exchanged  lithius  niobate", 
SPIE  1177.  Boston,  MA,  1989. 
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Figure  13.  Prooess  saquenoe  for  cncea^ing  proton-exciiaii^ 
waveguide  in  lithiue  nictate 
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Figure  13.  Scheeatic  of  an  electrode-less 
Mach-Sehnder  Interferometer 
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A  gradual  change  in  the  junction  region  could  be  obtained  by 
Bodifying  the  width  or  the  thickness  of  the  Ti  layer.  The  width  or 
the  thickness  could  be  tapered  over  a  length  of  a  few  hundred 
microns  at  the  junction  regions.  It  is  more  practical  to  taper  the 
width  than  the  thickness  of  the  Ti  layer. 

For  example,  if  the  desired  width  of  the  Ti  layer  is  94m,  the  width 
of  the  Ti  layer  could  increase  from  zero  (in  practice,  less  than 
0.25^m)  to  the  desired  width  of  9|im.  This  can  be  easily 
accomplished  by  incorporating  the  design  in  the  photomask  used  to 
define  the  Ti  pattern  on  the  suJt^trate. 

Since  the  domain  inversion  effect  depends  on  the  concentration  of 
titanium,  gradual  domain  inversion  would  also  occur  in  the 
transition  region.  A  sketch  of  such  a  device  that  has  a  tapered  Ti 
width  is  shown  in  Figure  14.  If  each  of  the  transition  regions 
required  0.5  mm,  and  the  total  length  of  the  Ti  region  is  10  mm, 
100%  domain  inversion  takes  place  only  over  9  mm  length  of  the  Ti 
diffused  region.  If  we  assume  that  the  transition  region 
contributes  to  only  50%  domain  inversion,  a  decrease  of  5%  in  total 
phase  shift  occurs  compared  to  a  non-tapered  design.  This  could  be 
accommodated  by  appropriately  increasing  the  straight  segment  to  be 
equal  to  9.5  mm. 
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Lithium  Niobito  Sutwtralt 


Figure  14.  Scheaatic  of  a  tapered  Ti  regirni  to 
overcoae  the  effects  of  abrupt  steps  in  the 
lower  leg  of  the  interferometer 
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An  electrode-less  Mach-Zehnder  interferometer  for  sensing 
electromagnetic  fields  is  shown  in  Figure  14.  The  Mach-Zehnder 
interferometer  is  fabricated  in  an  electro-optic  lithium  niobate 
crystal  material.  One  leg  of  the  interferometer  consists  of  2m 
ion-exchange  waveguide  within  a  Ti  diffused  domain  inverted  region. 
The  other  leg  of  the  interferometer  is  made  of  an  ion-exchanged 
waveguide.  There  is,  therefore,  a  built-in  asymmetry  in  this 
device  leading  to  a  phase  bias  in  the  absence  of  any  external 
electric  field. 

Mhen  a  Z-directed  electric  field  is  applied  to  the  device,  it 
changes  the  refractive  index  of  the  material  through  the  linear 
electro-optic  effect.  The  light  bean  experiences  a  phase  shift  as 
it  propagates  through  the  waveguide.  Since  one  leg  of  the 
interferometer  is  oppositely  poled,  it  experiences  an  equal  but 
opposite  phase  shift  to  that  of  the  other  leg  of  the 
interferometer.  This  device  thus  operates  in  the  sane  fashion  as 
an  electroded  device  described  earlier. 


The  quadrature  phase  bias  can  be  implemented  in  a  Mach-Zehnder 
using  several  approaches.  In  the  first  approach,  a  rc  bias  is 
applied  to  one  leg  of  the  interferometer.  The  DC  voltage  creates 
a  phase  difference  through  the  linear  electro-optic  effect.  Since 
this  requires  metal  electrodes  to  be  placed  on  the  crystal,  it  is 
not  applicable  in  the  case  of  a  non-metallic  electric  field  sensor. 

The  second  approach  is  based  on  creating  an  optical  path  length 
difference  in  the  device  design  and  fabrication  process.  The 
optical  path  length  is  a  product  of  the  physical  length  and  the 
effective  index  of  the  guided  mode  in  the  waveguide.  Hence,  a 
difference  in  the  physical  length  or  the  effective  index,  n^,, 
between  the  two  legs  of  the  interferometer  could  be  used  to  create 
a  net  phase  shift  of  r/2. 

In  the  case  where  the  two  legs  of  the  interferometer  have  lengths 
Li  and  L,  and  the  effective  index  is  the  same,  the  net  phase 
difference  in  the  absence  of  any  electric  field  is 

^  |£.  -  I^\  (17) 


or 

\A^,\  =  (18) 
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For  a  i^sa  diffarance  of  r/2,  the  physical  path  langth  diffaranca 
raquirad  is 

AL  .  A  ^  (19) 


For  X  >  1.3  iui  and  «  2.15» 

AL  >  0.15  liB. 

Such  a  swill  path  diffaranca  is  difficult  to  fabricata  using 
current  technology  of  lithography  and  processing.  Fortunately,  the 
periodic  nature  of  the  cosine  squared  function  associated  with  the 
MZ  transfer  curve  can  be  used  to  ispleeent  a  idiase  difference  that 
is  an  odd  aultiple  of  r/2  or  an  odd  aultiple  of  AL,  such  that 

AL  ■  X  0.15  \XM. 

The  fabrication  process  is,  therefore,  relaxed  by  choosing  a  path 
length  difference  of  several  nicrons. 

In  the  case  iirtiere  the  w/2  phase  bias  is  obtained  by  creating  two 
different  waveguide  effective  indices,  the  device  is  designed  such 
that  only  the  two  straight  sections  have  a  differing  n^.  The 
input  and  output  Y-branches  have  identical  effective  indices  so 
that  they  do  not  contribute  to  any  jdiase  bias. 

If  the  length  of  the  straight  section  is  L,  and  n^m  and  n^,  are 
their  effective  node  indices,  the  net  phase  bias  is 


lAAcI 


(20) 


|A^e|  *  ^^0tf 


(21) 


A  r/2  phase  bias  is  obtained  when 


(22) 


The  branching  angle  at  the  input  and  output  of  a  Nach>Zehnder 
inter feroseter  is  snail,  on  the  order  of  one  degree.  It  is,  thus, 
nore  convenient  to  keep  the  physical  size  of  this  interferoneter 
balanced.  To  do  this  and  still  acihieve  kmilt-in  quadrature  bias, 
we  will  first  try  to  achieve  the  quadrature  phase  shift  by  choosing 
the  length  of  the  ams  appropriately.  The  presence  of  Ti  in  one 
am  will  lead  to  an  inherent  phase  difference  between  the  two  ams. 
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If  it  is  possible  to  choose  the  am  length  so  as  to  give  a  phase 
difference  of  an  odd  multiple  of  w/2,  quadrature  operation  may  be 
achieved  in  this  way. 

The  quadrature  bias  condition  for  the  electrode-less  device  under 
consideration  could  be  obtained  using  the  optical  path  length 
difference  based  on  the  difference  in  the  effective  indices  of  the 
two  legs  of  the  interferometer.  The  electric  field  sensor  concept 
being  pursued  involves  a  Mach-Zehnder  interferometer  fabricated  on 
a  LiNbO,  substrate  in  which  one  am  consists  of  an  optical  channel 
waveguide  fomed  by  proton-exchange  and  one  am  fomed  by  proton- 
exchange  plus  in-diffusion  of  titanuim.  Because  both  Ti  diffusion 
and  proton-exchange  each  cause  an  increase  in  refractive  index,  the 
interferometer  is  unbalanced.  The  object  of  this  design  is  to 
allow  the  interferometer  to  operate  with  a  quadrature  phase  shift 
between  the  two  arms. 

The  cross  section  of  the  upper  leg  shown  in  Figure  15  consists  of 
a  proton-exchange  waveguide  in  a  LiNbO,  substrate.  This  is 
illustrated  in  Figure  16  as  a  three  layer  model  consisting  of  a 
lithium  niobate  substrate,  proton-exchange  waveguide  and 
superstrata  (air).  The  proton-exchange  increases  only  the 
extraordinary  refractive  index  of  the  crystal  and  does  not  affect 
the  ordinary  refractive  index.  Hence,  the  ordinary  refractive 
index  of  the  proton-exchange  waveguide  region  is  the  aama  as  that 
of  the  substrate. 

The  cross  section  of  the  lower  leg  shown  in  Figure  17  consists  of 
a  proton-exchange  waveguide  within  the  domain  inverted  region  of 
the  lithium  niobate  substrate.  This  is  illustrated  in  Figure  18 
as  a  three  layer  model  consisting  of  a  Ti  diffused  lithium  niobate 
substrate  region,  proton-exchange  waveguide  and  superstrata  (air). 
As  explained  earlier,  the  proton-exchange  process  does  not  affect 
the  ordinary  index  of  the  waveguide.  Hence  the  ordinary  refractive 
index  of  the  waveguide  region  is  the  same  as  that  of  the  Ti 
diffused  region. 

We  shall  now  analyze  the  difference  in  effective  refractive  index 
between  the  two  arms  of  the  interferometer  described  cUbove.  He 
begin  our  analysis  by  considering  a  three  layer  planar  structure  in 
which  the  three  layers  are  the  substrate,  the  waveguide  region  and 
air  having  refractive  indices  n,,  n,  and  n,,  respectively.  We  assume 
no  electromagnetic  field  variation  in  the  y  direction  and  that 
propagation  in  the  z  direction  is  described  by 

(23) 

where  B  is  the  mode  propagation  constant. 

The  electromagnetic  wave  equation  may  then  be  solved  for  the 
transverse  electric  (TE)  and  transverse  magnetic  (TM)  cases.  The 
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Proton-Exchange  Waveguide 
Subetrate 


Figure  15.  Croee  section  of  the  proton-eiffduinge  waveguide 
in  the  uppeac  leg  of  the  Macdi^Sehnder  interferoneter 
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Figure  16. 


Three  layer  eodtel  of  the  PE  waveguide 
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Figure  17.  Cross  sectim  of  the  proton-exchange  imveguide 
in  the  Ti  diffused  doauiin  inverted  regicni  in  the  lower  leg  of 
the  Macdi-Zehnder  interf eroaeter 
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Figure  18.  niree  layer  sodel  of  the  PB/Ti  waveguide 
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Figure  19.  Four  layer  sodel  of  the  PB/Ti  waveguide 
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(24) 


effective  refractive  index  is  defined  as 


n 


mff  • 


1 

K 


where  x  -  2w/X  and  X  is  the  free  space  wavelength. 

We  consider  solutions  of  the  rave  equation  for  n,  <  n^,  <  n,, 
assuaing  n^  >  n^.  Iiqposing  boundary  conditions  and  radiati<Mi 
conditions  at  x  *  ±oo  reguires  that  x  for  each  node  sust  satisfy 

tan  ht  «  ^  (25) 

for  TE  and  TM  aodes,  where 


g  «  (P*  - 

(26) 

h  - 

(27) 

P  «  (P»  - 

(28) 

and 


0^-1 


( TB  nodes) 


(29) 


o 


i 


(JM  nodes) 


(30) 


-  OjOjpgr 

(31) 

OjJi  (o^g  +  O3P) 

(32) 

We  have  used  the  above  equations  to  evaluate  the  effective 
refractive  index,  n^,  for  the  two  cases  of  (1)  a  proton^exchanged 
waveguide  and  (2)  a  proton~exchanged  waveguide  in  idiich  Ti  has  also 
been  diffused.  A  Fortran  prograa  idiich  uses  double  precision  has 
been  used  to  sake  this  calculation.  We  have  considered  single  node 


DAAL02-90-C-0078 
Optical 


45 

for  Sensing 


8RIOO 

Fields 


waveguides  and  the  TM  polarization  for  the  planar  waveguide,  ffe 
have  Bade  realistic  assuaptions  that  the  presence  of  Ti  causes  the 
refractive  index  to  increase  by  0.001  and  that  the  proton-exchange 
process  leads  to  a  refractive  index  increase  of  0.003  for  the  TM 
node.  In  this  calculation  we  assune  that  the  Ti  diffusion  is 
sufficiently  deep  so  that  the  diffused  region  can  act  as  the 
substrate.  Later,  we  detemine  the  depth  for  which  this  assumption 
is  valid.  The  parameters  used  and  the  results  of  the  calculation 
are  listed  in  Table  III.  The  change  in  effective  index  due  to  the 
presence  of  Ti  is  0.001001  or  just  0.001,  the  same  value  as  the 
actual  refractive  index  change  we  assumed.  This  sane  result  has 
been  obtained  for  a  range  of  parameters,  that  is,  the  change  in 
effective  index  is  equal  to  the  change  in  actual  index. 

The  change  in  phase  A*  due  to  the  above  change  in  effective  index 
An«f(  for  a  Mach-Zehnder  interferometer  am  length  L  is  given  as 

(33) 


Quadrature  operation  implies  that  this  phase  shift  is 

=  in-~  /B  =  1,3.5, . 

2 


(34) 


This  can  be  accomplished  for  lengths  L  given  as 


L 


»  m 


X 

4An.„ 


(35) 


Using  An«ct  »  0.001  from  above  and  X  -  1.3  un,  we  find 

L  »  m  X  324.7  micron  m  »  1,3,5,.... 

This  is  a  small  enough  value  for  n  »  1  that  achieving  an  odd 
multiple  of  it  is  realistic.  The  only  question  is  the  accuracy  of 
the  above  calculation. 

We  can  improve  our  model  by  using  a  4-layer  planar  fomulation  on 
the  interferometer  am  which  has  undergone  both  proton-exchange  and 
Ti  diffusion.  By  considering  an  extra  layer  we  can  examine  the 
effect  of  varying  the  depth  of  the  Ti  diffusion  has  on  n^^^.  For 
the  3-layer  calculation  above,  we  assumed  that  this  depth  was 
sufficiently  deep  for  the  Ti  diffused  layer  to  act  as  the 
substrate.  The  waveguide  parameters  for  the  4-layer  fomulation 
are  shown  in  Figure  19. 

The  thickness  of  the  waveguide  is  t  and  the  depth  of  the  Ti 
diffused  region  is  d,  with  the  difference  6  given  as 
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3>Layer  Calculation  of  Effective  Indices 


Case  1 

Proton  Exchange  (PE) 

Case  2 

PE  Ti-Indiffusion 

atr 

1.0 

1.0 

^guUe 

2.153 

2.154 

^substrata 

2.150 

2.151 

t  (im) 

7.75 

7.75 

X(iim) 

1.3 

1.3 

^Single  mode 

8.56 

8.56 

a  planar 

2.151960 

2.152961 
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6  *  d  -  t 


(36) 


Th«  •ffttctive  indax,  B,  for  the  above  4-layer  structure  is  given  as 
(o^g  ♦  Oqr)  (|ft  -  Ytan  bt)  *  (OgZ  -  o^g)  x 

♦OiOsPg)  tanht  ♦  a^io^g  -  OjP)]  »  0  (36) 

with  the  BOde  paraaeters  defined  in  Equations  (25)  -  (31)  and 

y  -  (f>  -  It>n,»)^ 


We  have  calculated  the  effective  refractive  index  n«,,(6)  as  a 
function  of  S,  defined  in  Eg.  (35).  If  6  »  0,  then  the  Ti  depth  is 
the  saM  as  the  proton-exchange  depth.  As  6  increases,  n«M(6)  also 
increases.  For  convenience,  we  plot  this  increase  relative  to  the 
value  of  n«ff  at  6  0  in  Figure  20.  Note  that  this  difference 
increases  with  S  initially,  but  that  it  saturates  once  4  «  2  tm. 
Figure  12  provides  t%ro  iaportant  pieces  of  inforaation:  (1)  the 
value  of  S  for  which  the  difference  saturates,  and  (2)  the  value  of 
this  difference  at  saturation.  Although  the  actual  index 
difference  at  saturation  is  small,  »  6  x  10“*,  the  corresponding 
Idtiase  shift  for  a  10  n  long  interaction  length  is 

(10  X  10’X«-0  X  10'*)  •  0.9  *  (38) 

which,  of  course,  is  significant.  In  the  above  equation,  6n^ 
represents  the  change  in  effective  index,  plotted  in  Figure  12,  in 
the  lower  leg  of  the  interferoaeter  due  to  finite  6.  Thus,  if  6 
varies  between  zero  and  the  value  for  idiich  the  change  in  Figure  20 
saturates,  »  2  |ui  in  this  example,  it  will  have  a  significant 
effect  on  determining  the  quadrature  condition.  Because  the  exact 
value  of  6  may  be  hard  to  determine  or  hard  to  control,  the  best 
procedxire  may  be  to  insure  that  the  Ti  diffusion  depth  is 
sufficiently  large  so  that  the  change  in  effective  index  6n^,  shown 
in  Figure  20  has  saturated.  For  this  example,  the  Ti  diffusion 
depth  would  need  to  be  at  least  2  iim  greater  than  the  proton- 
exchange  depth. 

We  can  refine  the  calculations  for  the  example  considered  above  by 
including  the  effect  of  the  channel  waveguide.  The  simplest  %my  to 
do  this  is  to  use  the  effective  index  method.  To  malce  this 
calculation,  it  is  assumed  that  the  Ti  diffused  region  is  much 
wider  than  the  proton-exchange  region.  The  results  of  calculating 
the  channel  %raveguide  effective  index  are  presented  below  in 
Table  IV.  Note  that  although  the  effective  indices  are  different 
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Efifecttve  Index  Calculation  of  Channel 


Waveguide  Effective  Indices 


Case  1 

Proton  Exchange  (PE) 

Case  2 

PE  +  Ti'Indiffusion 

e  planar 

2.151960 

2.152961 

dad 

2.150 

2.151 

w  (fim) 

7.0 

7.0 

X  (tun) 

1.3 

1.3 

w  ,  ,  _ _  (tun) 

dngto  mode  ' 

7.08 

7.08 

n 

e 

2.151259 

2.152260 
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than  those  calculated  in  Table  III,  the  difference  in  effective 
indices  is  the  saee.  This  implies  that  the  value  of  length 
calculated  from  Eq.  (34)  needed  to  achieve  quadrature  offset  would 
be  the  same  as  determined  earlier,  that  is 

L  »  m  X  324.7  iim  m  >■  1,3,5,.... 

This  value  would  be  modified  by  4-layer  considerations,  as 
discussed  earlier. 

^ese  calculations  show  that  it  may  be  possible  to  choose  the 
quadrature  point  of  operation  by  appropriately  choosing  length,  lie 
could  improve  this  possibility  by  repeating  the  above  calculations 
for  actual  graded  index  profiles,  assuming  we  could  predict  the 
proton-exchange  and  Ti  diffused  index  profiles  accurately.  Because 
small  effects,  such  as  the  4-layer  effect  discussed  earlier,  can 
significantly  alter  the  quadrature  point,  it  may  still  be  difficult 
in  practice  to  theoretically  predict  the  quadrature  point,  even  if 
elaborate  graded  index  approaches  were  used.  Nhat  is  needed  is 
some  way  to  adjust  the  quadrature  point  once  a  device  has  been 
fabricated. 

Tuning  the  source  wavelength  may  be  one  such  way  to  adjust  the 
phase  difference  between  the  two  arms  of  the  interferometer.  We 
consider  such  a  possibility  for  the  parameters  given  in  Table  III 
and  for  an  interaction  length  of  10  mm.  We  see  in  Figure  21  that 
the  difference  in  the  effective  indices  of  the  two  arms  of  the 
interferometer  is  nearly  constant  with  wavelength.  Although  the 
effective  index  in  each  arm  varies  with  wavelength,  the  variations 
in  each  arm  are  similar  so  that  the  difference  remains  constant. 
However,  we  see  in  Figure  13  that  significant  phase  change  occurs. 
This  phase  change  is  due  to  the  2w/k  term  in  Eq.  (32).  For  the 
range  of  wavelengths  considered  in  Figure  13  the  i^se  change  is 
more  than  sufficient  (±r/2)  to  choose  a  value  of  %ravelength  for 
quadrature  operation.  The  calculation  in  Figure  13  is  for  L  > 
10  am;  for  longer  lengths  the  range  of  wavelengths  would  be 
smaller.  Specifically,  to  achieve  a  phase  variation  of  ±ir/2,  the 
range  of  wavelength  over  which  the  source  must  be  tuned  is 

measured  in  fractions  of  micron.  For  example,  if  L  »  20  am, 
“  *001  and  A  »  1.3  mm,  then  the  range  of  wavelength  over  which 
the  source  needs  to  be  tuned  is  0.034  |ia  to  allow  variation  of  the 
fhase  by  ±w/2,  thus  assuring  quadrature  operation  can  be  achieved. 

The  discussion  here  has  been  directed  at  establishing  a  basis  for 
designing  the  asymmetric  Mach-Zehnder  interferometer  so  that 
quadrature  phase  difference  between  light  emerging  from  the  two 
arms  can  be  achieved.  Numbers  used  here  are  representative  of  the 
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Figortt  21.  Plot  of  the  cdiango  in  offoctive  indox  botwoon  tlM  two 
urmm  ot  the  interf eroMter  and  the  interf  eroaeter  phase  change  ae 
a  function  of  wavelength  over  a  United  range 
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technologies,  and  results  of  the  calculations  presented  here 
provide  a  basis  for  actual  design.  Planar  3-layer  and  4-layer  and 
channel  optical  waveguide  dispersion  calculations  have  been 
performed.  These  calculations  predict  a  length  for  quadrature 
operation  which  could  be  achieved  in  practice,  but  these  calculated 
lengths  are  sensitive  to  the  refractive  index  variations  associated 
with  the  proton-exchange  and  Ti  diffusion  technologies.  Although 
more  sophisticated  models  could  be  utilized,  it  is  expected  that 
their  use  would  still  result  in  the  same  problem.  We  thus  have 
investigated  the  possibility  of  using  wavelength  tuning  to  adjust 
the  phase  difference  for  quadrature  operation.  The  range  of 
wavelength  tuning  required  is  sho%m  to  be  realistic. 

Another  approach  is  to  modify  the  depth  of  the  proton-exchange 
waveguide  within  the  Ti  diffused  region.  The  effective  index  is  a 
sensitive  function  of  6  as  shown  in  Figure  20.  We  could  use  this 
property  of  the  waveguide  to  tailor  the  bias  condition  using  post¬ 
processing  techniques.  Since  the  proton-exchange  waveguide  could 
be  driven  deeper  at  a  relatively  low  temperature  of  200  to  400  C, 
a  strip  heater  could  be  used  to  heat  the  device  to  change  the  bias 
condition.  The  device  output  could  be  monitored  during  the  heating 
process  and  the  heater  turned  off  once  the  quadrature  bias 
condition  is  reached.  This  could  be  implemented  during  the  final 
testing  and  packaging  stages  of  the  sensor. 
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LiNbO,  is  transparent  from  the  visible  (0.45  ^m)  to  the  mid- 
infrared  (4  ^m) .  In  principle,  the  device  could  operate  anywhere 
in  the  transmission  spectrum  of  the  crystal.  The  device  under 
consideration  is  based  on  integrated  optical  technology  which 
requires  single  mode  optical  fibers  for  optical  input  and  output. 
Single  mode  optical  fibers  are  commercially  available  for  operation 
in  the  near  infrared  wavelength  range,  we  are  limited  to  operating 
this  device  within  this  wavelength  band.  The  wavelength  range  is 
divided  into  two  bands  —  0.8  to  0.9  |im,  the  short  wavelength 
region,  and  1.2  to  1.6  ^m,  the  long  wavelength  region.  Excellent, 
high  quality,  solid  state  laser  diodes  are  also  commercially 
availzd>le  at  these  wavelengths. 

Because,  as  will  later  be  seen,  several  milliwatts  of  optical  power 
may  have  to  be  used  to  achieve  the  desired  dyneuaic  range,  a 
wavelength  of  1300  nm  instead  of  850  nm  is  preferred,  so  as  to 
reduce  or  eliminate  optical  damage  effects  in  LiNbO,. 

The  selection  of  wavelength  of  operation  is  based  on  the  following 
considerations : 


Device  sensitivity 

Optical  power  damage  threshold 

Device  fabrication  and  packaging 


DEVICE  SENSITIVITY 

The  total  phase  shift  in  a  Nach-Zehnder  interferometer  is  given  by 

|A*.|  =  2  -J  nj  r„  1,  (40) 


The  phase  shift  is  inversely  proportional  to  the  wavelength  and 
directly  proportional  to  the  cube  of  the  refractive  index.  In 
general,  the  variation  of  refractive  index  with  wavelength  is  very 
small  and  we  can  thus  conclude  that  the  total  phase  shift  decreases 
as  wavelength  increases.  The  signal  to  noise  ratio  (SNR),  hence 
the  device  sensitivity,  is  proportional  to  the  square  of  the  phase 
shift. 


The  improvement  in  device  sensitivity  for  a  device  operating  at  0.8 
fim  and  1.3  ^m  is  given  by 


(  A«|>o.a  V 

1  j 

=  (^L 

M.. 

(41) 
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For  X  «  0.8  iui,  n«  »  2.17,  and  X  «  1.3  lua,  n,  «  2.15, 

»  2.8a4.5dB  (42) 


We  s«e  that  a  device  operating  at  0.8  im  %K>uld  give  an  iqprov«Mnt 
in  device  sensitivity  of  4.5  dB  based  on  the  above  calculations. 
This  is  based  on  the  assusption  that  the  device  is  capable  of 
carrying  the  saae  anount  of  optical  power  at  both  these 
wavelengths.  This  is  not  true  in  LillbO,  due  to  optical  power 
dasage  constraints. 


OPTICAL  POWBR  DMaCE  imKgHQTJl 

Optical  dasage  is  also  known  as  {diotorefractive  effect  trhere  the 
refractive  index  of  the  Material  is  dependent  on  the  intensity  of 
the  optical  beaa.  The  optical  power  density  at  which  the  Material 
exhibits  this  effect  is  knoim  as  the  threshold  of  optical  dasage. 

Optical  dasage  occurs  in  a  LiNbO,  crystal  because  of  the  presence 
of  iMpurities.  The  crystal  received  froM  the  Manufacturer  has  less 
than  a  few  pps  of  iron  and  other  rare  earth  ions.  Additional 
inpurities  are  introduced  into  the  crystal  during  the  device 
fabrication  process.  For  exaaple,  the  Ti  aetal  used  in  diffusion 
carries  with  it  additional  iron  inpurities  such  that  a  diffused 
region  of  the  substrate  could  have  tens  of  ppn  of  iron  present. 

The  effect  of  an  optical  bean  in  the  crystal  is  to  ionize 
inpurities  such  as  iron.  The  optical  field  in  the  crystal 
separates  the  free  carriers  fron  the  inpurities  and  traps  these 
carriers  in  unilluninated  regions  resulting  in  a  local  electric 
field.  This  local  field  changes  the  refractive  index  in  that 
region  through  the  linear  electro-optic  effect  and  alters  the 
waveguide  characteristics  of  the  crystal.  The  overall  effect  is 
that  the  optical  wave  initially  confined  in  the  waveguiding  region 
is  radiated  into  the  substrate.  The  extent  of  optical  daaage  is 
dependent  on  the  power  density  and  the  wavelength  of  the  optical 
signal  in  the  waveguide.  Optical  power  dasage  is  significant  at 
lower  wavelengths  such  as  0.8  lua  coMpared  to  1.3  |ur. 

The  optical  waveguide  in  LiNbO,  is  capedsle  of  handling  10  to  100  \M 
of  optical  power  at  the  wavelength  of  0.8  ixn  and  greater  than  10  mH 
of  optical  power  at  1.3  ^m  wavelength.  To  obtain  the  required 
device  sensitivity,  we  would  require  to  propagate  several 
Milli«ratts  of  optical  power  in  each  waveguide  channel.  Based  on 
the  optical  dasage  considerations,  we  conclude  that  the  wavelength 
of  operation  should  be  1.3  The  Mach-Zehnder  device  optisized 

to  operate  at  1.3  im  could  also  be  operated  at  1.55  |ib  wavelength. 
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■niis  would  provide  increased  optical  power  handling  capability  in 
the  waveguide  but  a  decreased  swtsitivity  of  1.5  dB. 


pgvici  rmicmoM  MP  PMaaisinG 

The  waveguide  linewidth  at  1.3  lui  and  0.8  lis  wavelengths  are  7  tm 
and  4  im  respectively.  The  larger  linewidth  saXes  the  task  of 
device  fabrication  easier  at  1.3  and  also  increases  the 
production  yield.  Similarly,  the  mode  size  of  the  optical 
waveguide  and  the  fiber  core  size  are  larger  at  1.3  im  wavelength. 
The  alignment  accuracy  between  an  optical  fiber  and  the  device  at 
0.8  |Ui  is  in  the  submicron  range.  A  id.salignsent  of  0.5  luz  reduces 
the  optical  coupling  efficiency  betiraen  the  fiber  and  the  traveguide 
by  3  dB.  The  waveguide  device  at  1.3  \m  is  extremely  tolerant  to 
misalignment  of  at  least  3  tim  before  exhibiting  a  3  dB  optical 
coupling  loss. 

Based  on  these  considerations,  the  wavelength  of  operation  is 
chosen  to  be  1.3  4m. 
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This  section  provides  an  analysis  of  the  proposed  optical  sensor 
systeo.  The  design  of  the  electric  field  sensor  is  influenced  by 
the  interrelationship  of  several  systen  paraneters.  The  carrier- 
to-noise-ratio  (CNR)  for  a  given  received  optical  power  Bust  first 
be  determined.  From  this,  the  signal-to-noise-ratio  (SNR)  for  a 
given  peak  phase  modulation  index  is  derived.  This  peak  phase 
modulation  index  limits  the  dynamic  range  and  is  based  on  linearity 
considerations.  The  minimum  detectable  electric  field  is 
calculated  based  on  SNR*b1.  The  maximum  detectable  electric  field 
is  then  calculated  using  the  device  geometry  and  the  crystal 
properties.  The  minimum  and  maximum  electx'ic  field  yield  the 
dynamic  range.  The  results  of  these  calculations  are  summarized  in 
a  tabular  form  at  the  end  of  this  section. 


The  carrier-to-noise-ratio  (CNR)  in  a  single-ended  homodyne 
detection  optical  system  is  given  by  the  following  relationship: 


neP/i) 


CNR  = 


,(3g. .  .  iSt. 


(43) 


where  P.  -  received  optical  power  in  each  beam  as  if  acting 

alone, 

r\  =  quantum  efficiency, 
h  =  Planck's  constant  (6.63  x  10"”  J.s), 

=  optical  carrier  frequency  (2.31  x  10”  Hz  for  a 
1300  nm  source), 

=  dark  current, 

M  =  avalanche  gain  («  1  for  a  PIN  photodiode) , 

X  =  excess  noise  factor  (=0  for  a  PIN  photodiode), 
k  =  Boltzmann's  constant  (1.38  x  10"**  J/K), 

F  =  noise  figure  of  eusplifier, 

Rl  ==  effective  load  resistance, 

B.  =  noise  equivalent  bandwidth  (1.57  GHz  for  B.  equal  to 
a  1  GHz  signal  system  bandwidth). 

In  a  system  where  a  single  pole  low  pass  filter  defines  the 
overall  bandwidth,  the  noise  equivalent  bandwidth  is  x/2  times  the 
-3dB  system  bandwidth. 


Substituting  the  values  of  P. 
the  above  equation. 


1  mW,  1)  =  0.67  for  a  PIN  diode  into 


CNR  =  63.9  dB. 
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In  the  Equation  (43),  the  CNR  increases  when  the  optical  power,  P., 
increases  or  when  the  bandwidth,  B.,  decreases.  Increasing  the 
load  resistor,  R,.,  will  decrease  the  theraal  noise,  but  it  will 
lower  the  frequency  response.  Moreover,  it  will  render  the 
detector  unable  to  detect  high  frequency  signals.  The  avalanche 
multiplication  factor  M  tends  to  increase  the  CNR,  particularly 
when  the  optical  power  is  very  low  and  the  thermal  noise  is 
dominant.  The  multiplication  factor  (M)  for  a  PIN  diode  is  1, 
while  that  for  an  APD  is  typically  between  10  and  100. 

To  maximize  the  sensitivity  and  dynamic  range,  the  optical  receiver 
is  designed  to  be  shot  noise  limited.  The  shot  noise  or  the  kT 
noise  at  300  K  is  -174  dBm/Hz.  If  the  effective  photodetector  load 
was  50  ohms,  and  the  following  eunplifiers  were  essentially 
noiseless,  i.e.,  noise  figure  F  =  1,  then  the  power  reaching  the 
photodetector  would  have  to  be  sufficient  to  ensure  that  the  shot 
noise  generated  in  the  load  would  be  >  -174  dBm/Hz.  Anything  less 
than  that  would  limit  the  CNR  to  a  value  less  than  optimum.  By 
using  a  transimpedance  amplifier  with  a  PIN  photodetector,  the  load 
can  be  increased  to  2dx>ut  IK  ohm  without  compromising  bandwidth. 
Thus,  it  reduces  the  power  needed  to  reach  the  quantum  noise  limit. 

The  required  dynamic  range  is  40  dB,  corresponding  to  phase  and 
electrical  field  ratios  of  100:1.  This  would  be  adequate  for 
sufficiently  large  phase  modulation.  In  most  cases,  the  peak  phase 
modulation  A#  may  have  to  be  considerably  less  than  unity  in  order 
to  ensure  a  low  level  of  distortion  from  the  phase  detector.  The 
SNR  could  fall  below  40  dB  in  such  situations. 

The  optical  source  could  also  degrade  the  SNR  due  to  excess 
intensity  noise.  This  general  broadband  excess  noise  level  is 
commonly  expressed  as  the  Relative  Intensity  Noise  (RIN) .  The  RIN 
generally  has  a  fairly  flat  spectrum,  except  at  (a)  frequencies 
determined  by  external  cavity  resonances  due  to  reflections  back 
into  the  lase  diode,  (b)  near  the  laser  relaxation  frequency  and 
(c)  near  zero  hertz.  Since  excess  noise  normally  varies 
substantially  with  optical  feedback,  the  use  of  an  optical  isolator 
is  recommended.  Optical  isolators  are  more  compact  and  less 
expensive  at  1300  nm  than  at  800  nm  wavelength. 

As  a  rule  of  thumb,  in  a  direct  detection  system  a  shot  noise 
limited  1  mW  near  infrared  optical  system  will  produce  a  CNR  of 
approximately  150  dBc/Hz  (dBs  with  respect  to  the  carrier  in  a  1  Hz 
bandwidth).  This  assumes  that  there  is  no  excess  noise.  If  the 
RIN  value  for  this  device  is  -140  dBc/Hz,  so  that  the  CNR  saturates 
at  high  power  levels  to  a  value  of  140  dB/Hz,  then  the  excess 
intensity  noise  is  about  10  dB.  The  dynamic  range  of  a  system 
using  this  device  will,  therefore,  be  degraded  by  about  10  dB. 

In  the  following  pages  are  included  Tables  V  and  VI  generated  from 
the  Optoelectronic  Spreadsheet  Template  CNR.NRl  developed  by 
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Table  V  —  CMR  for  slnale-ended  detection 
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PilMrdyiM  Optoelactronics.  The  highlighted  items  are  the  more 
im|>ortant  parameters  that  affect  the  cam.  Table  V  gives  the 
analysis  for  a  single-ended  detection  system.  For  a  negligible 
excess  noise  of  -200  dBc/Hz,  as  shown  in  Range__A,  the  CMR  is 
63.9  dB  for  a  received  optical  power  of  1  nw. 

The  CNR  drops  to  50.8  dB  for  an  excess  noise  of  -140  dBc/Hz.  The 
corresponding  grafA,  Figure  22,  shows  this  reduction  in  CNR  and 
the  dynamic  range  in  a  more  dramatic  %ray.  This  reduced  CNR  may 
present  a  major  barrier  in  obtaining  the  desired  dynamic  range. 
The  effect  of  intensity  noise  could  be  decreased  using 
com>lementary  idiotodetection.  To  do  this,  the  Mach-Zebmler 
interferometer  %iould  require  t%io  output  ports  instead  of  one.  This 
is  accomplished  by  the  addition  of  a  directional  coupler  at  the 
output  of  the  nach-Zehnder  as  shown  in  Figure  23. 

The  signals  applied  to  dual  i^iotodetectors  are  coiqplementary.  As 
the  phase  in  one  branch  of  the  interferc»eter  varies,  the  intensity 
on  one  photodetector  increases  tdiile  on  the  other  photodetector  it 
decreases.  Both  detectors  have  the  same  received  intensity  noise 
if  the  interferometer  is  in  phase  quadrature.  Thus,  %rhile  the 
difference  signals  between  the  irfiotodiodes  due  to  the  i^ase-induced 
intensity  changes  add,  since  they  are  eqpial  and  opposite,  the 
common  mode  intensity  noise  cancels.  Typically,  with  careful 
design,  up  to  30  or  more  dBs  of  intensity  noise  suppression  can  be 
obtained.  This  dual-detection  receiver  also  has  the  advantage  that 
the  half  of  the  optical  power  that  was  previously  wasted  in  the 
substrate  is  now  available  at  the  additional  photodetector.  Thus, 
while  the  electrical  signal  power  detected  by  the  front-end  is  four 
times  larger,  the  shot  noise  is  twice  as  large,  leading  to  a  net 
improvement  in  CNR  of  3  dB. 


Table  VI  gives  this  analysis  for  a  complementary  detection  system 
providing  30  dB  of  intensity  noise  suppression.  For  negligible 
excess  noise,  the  CNR  now  increases  by  3  dB  to  66.9  dB.  Even  when 
the  RIN  is  significant,  the  CNR  is  only  slightly  reduced  by  0.1  dB. 
This  is  practically  identical  to  the  results  obtained  from  the  shot 
noise  limited  expression  shown  below. 


am 


2neP, 

bfcBa 


(44) 


This  indicates  that  the  optical  front-end  is  shot  noise  limited  for 
P.  «  1  mW.  The  accompanying  grairii.  Figure  24,  demonstrates  this 
more  clearly.  Single  mode  laser  diodes  are  likely  to  have  less 
intensity  noise  because  of  the  lack  of  mode  competition  effects. 
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CNR  IN  HOMODYNE  OPTICAL  SYSTEMS 

Single-Ended  Detection 


RIN«-200dBc/Hz  RIN--160dBc/Hz  RIN  > -ISOdBc/Hz  RIN  - -140  dBc/Hz 
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CNR  IN  HOMODYNE  OPTICAL  SYSTEMS 

Dual-Balanced  Detection 


RiN>-200dBc/Hz  RIN  - -leOdBo/Hz  RIN  - -150  dBc/Hz  RIN  - -140  dBcAlz 
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This  section  deals  with  the  harsonic  distortion  caused  by  exciting 
the  electric  field  sensor  two  frequencies  and  e,  of  equal 
amplitudes.  The  second  and  third  harmonics  and  the  intermodulation 
frequency  will  be  analyzed.  The  dynamic  range  is  limited  by  the 
intermodulation  term  because  it  produces  a  higher  level  of 
distortion  than  the  third  harmonic. 

The  input/output  relationship  for  the  electric  field  sensor  is 
given  by 

P  1 

[1  +  cos(  +  ^.(sin  o>jt  +  sin 


where  »  bias  point  of  the  modulator, 

«.  s  peak  phase  deviation  produced  due  to  the  field. 

This  expression  can  be  expanded  with  a  power  series  into  individual 
harmonic  and  intermodulation  terms.  The  resulting  modulation 
ratios  are^*: 

Modulation  ratio  at  the  fundamental  frequency,  e: 


(46) 


Modulation  ratio  at  the  second  harmonic  frequency,  2u: 


(47) 


Modulation  ratio  at  the  third  harmonic  frequency,  3»: 


**S.A.  Wilcox  and  M.D.  Drake,  "Practical  system  design 
considerations  for  wideband  fiber  optic  links  using  external 
modulators",  SPIE  £12,  pp.  234-239  (1988). 
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(6 

24 

(-1)  *  *‘♦•’1 
'  ’  16 

sin  3ut 

1  *  cos4 

2  32 

Interaodulation  ( 2(»,  -  ) : 


L 

24 

sin  (2«»2-o^)  t 

1  *  cos^ 

L  ^  3  (♦.)*] 

2  32  . 

(49) 


Th«s«  Bodulation  ratios  are  plotted  in  Figures  25  and  26  as  a 
function  of  aodulator  bias  point  for  a  peak  |8uise  deviation  of 
0.10  and  0.25  radian.  At  guadrature  bias,  the  second  harnonic 
disappears  and  the  dynamic  range  is  limited  only  by  the 
interiwdulation  distortion.  As  the  bias  point  deviates  from 
quadrature,  the  second  harmonic  term  rises  rapidly.  Referring  to 
the  graph  for  0.25  radian,  the  second  harmonic  rises  to  the  level 
of  intermodulation  for  a  bias  drift  of  7  degrees.  Beyond  this,  the 
second  harmonic  term  severely  degrades  the  dynamic  range.  For  a 
peak  deviation  of  0.1  radian,  the  second  harmonic  rises  to  the 
level  of  the  intermodulation  term  for  a  bias  drift  of  3  degrees. 
As  the  peak  deviation  increases,  there  is  less  demand  on  the 
stability  requirements  for  quadrature  bias. 

The  peak  phase  deviation  allowed  will  depend  on  the  criterion  for 
distortion.  The  criterion  selected  to  limit  the  peak  phase 
deviation  is  when  the  intermodulation  distortion  term  equals  the 
noise  floor  level.  Based  on  Eqs.  (46)  -  (49),  the  modulation 
ratios  of  the  fundamental  and  the  intermodulation  terms  are  plotted 
at  quadrature  bias  as  a  function  of  peak  phase  deviation. 
Figure  27  shows  that  the  noise  floors  for  optical  powers  of  1  and 
10  mW  are  horizontal  lines.  The  graph  indicates  that  the  maximum 
phase  deviations  allowed  are  0.172  and  0.115  radian  for  optical 
powers  of  l  and  10  mW,  respectively. 

If  small  bandwidth  such  as  10  kHz  is  used  to  analyze  the  spectrum 
of  the  electromagnetic  energy  received,  there  is  no  dynamic  range 
problem.  However,  to  "capture"  the  signature  of  a  complex  random 
electromagnetic  waveform  containing  frequency  components 
approrching  1  GHz,  a  detection  bandwidth  of  1  GHz  would  be 
required. 
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Modulation  Ratio  (dB) 


F±^ia3re  2  5 


Theoretical  Modulation  Ratio 


Peak  Deviation  =  0.1  rad 
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Modulation  Ratio  (dB) 


Fi^uxTe  2  6 


Theoretical  Modulation  Ratio 
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Peak  Deviation  =  0.25  rad 


Modutalion  Ratio  (d^ 


Plgure  27 


Interferometer  Response  -  Quadrature  Bias 
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The  relationship  among  signal-to-noise  ratio  (SNR),  phase 
modulation  index  (•«)  and  the  carrier-to>noise  ratio  (CNR)  is: 

SNR  =  cm  (50) 


All  phase  deviations  are  peak  values. 

For  P.  »  1  mW,  i.e.,  >  2  mW,  and  a  maximum  phase  deviation  of 

0.172  radian,  expressed  in  dB  format: 

SNR  s  CNR  -  18.3  dB 

a*  48.6  dB 

For  P.  »  10  mW,  i,e,  P^.  >  20  mW,  and  a  maximum  phase  deviation  of 
0.172  radian: 

SNR  »  CNR  -  21.8  dB 
=  55.6  dB 

For  P,  »  10  mW,  the  CNR  is  77.4  dB  for  a  1.57  GH2  noise  equivalent 
bandwidth.  The  minimum  detectable  phase  deviation  is  obtained  when 
the  SNR  =  0  dB.  For  a  1.57  GHz  noise  equivalent  bandwidth,  the 
minimum  detectable  peak  phase  deviation  is: 

f o  ^  =  1.9  X  10‘*  radian 

For  P.  -  1  mW,  and  a  1.57  GHz  noise  equivalent  bandwidth, 

^  =  6.4  X  lO"*  radian 

An  electric  field  40  dB  above  these  values  required  to  produce  a 
signal  at  this  noise  floor  level  would  produce  a  peak  phase 
deviation  several  dBs  below  the  maximum  phase  deviation  calculated 
above. 


iyrjy,  *  cm 


Analyses  of  dynamic  range  are  shown  in  Figures  28  and  29  for  a 
dual-detection  system.  The  SNR  penalty  produced  by  phase 
deviations  of  0.172  and  0.115  radian  for  1  and  10  mW  received  per 
channel  and  the  noise  floors  for  the  different  RIN  levels  are  shown 
as  horizontal  lines. 
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Figure  28 
Relative  Signal  Levels 
For  1  mW  Dual-Balanced  Detection 
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Figure  29 
Relative  Signal  Levels 
For  10  mW  Dual-Balanced  Detection 
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At  a  power  level  below  1  nW,  it  nay  be  difficult  to  maintain  shot 
noise  limited  detection  consistent  with  the  required  1  GHz  system 
bandwidth  and  with  the  use  of  PIN  photodetectors.  The  graphs 
produced  from  Tables  V  and  VI  show  that  below  P.  «  0  dBm  (1  mW) , 
the  optical  receiver  starts  to  become  kT  noise  limited.  A  good 
rule  of  thumb  for  a  high  quality  PIN  photodetector  is  that  for  a 
photodetector  load  of  Ik  ohm,  ed)out  1  mH  of  total  power  (2P.)  is 
required  to  achieve  the  threshold  of  quantum  noise  limited 
detection.  It  is  not  possible  to  increase  the  load  much  above  this 
value  without  compromising  the  bandwidth.  It  would  be  better  to 
avoid  using  APDs,  particularly  for  the  dual -detect! on  system.  At 
these  reduced  power  levels  it  probably  becomes  less  of  a  necessity 
to  consider  complementary  photodetection,  unless  the  excess  noise 
is  very  large. 

The  choice  of  single  or  dual-detection  depends  on  the  RIN  values 
for  commercial  1300  nm  lasers.  If  a  complementary  photodetection 
system  is  not  used,  a  laser  with  a  cut  off  frequency  considerably 
in  excess  of  1  GHz  will  be  required.  This  must  be  done  to  ensure 
that  the  low  frequency  wing  of  the  relaxation  peak,  which  can  be 
accompanied  by  more  than  80  dB/Hz  of  excess  noise,  does  not  raise 
the  noise  floor  near  the  1  GHz  system  -3  dB  frequency. 


SENSITIVITY  ANALYSIS 

There  are  two  main  differences  between  the  electroded  and  the 
electrode-less  device.  First,  the  dielectric  constant  of  lithium 
niobate,  is  greater  than  that  of  free  space  and,  thus,  tends  to 
"short  out"  the  electric  field.  Second,  the  conventional  dipole 
antenna  has  a  voltage  induced  on  it  which  is  then  applied  to 
closely  spaced  electrodes  on  a  substrate.  A  very  large  localized 
field  is  thus  produced  in  an  electroded  device.  Hence,  the 
electric  field  "seen"  by  the  electrode-less  Nach-Zehnder 
interferometer  will  be  substantially  less  than  for  the  device  with 
electrodes . 

This  poses  no  problem  for  detecting  Electro  Magnetic  Pulse  (EMP), 
which  mainly  covers  the  spectral  range  10  kHz  to  100  MHz  and  has 
peak  field  strengths  of  50  kv/m.  Indeed,  the  presence  of  discrete 
metallic  electrodes  could  cause  electrical  breakdown  of  the  device. 
Without  electrodes,  the  electric  field  may  be  allowed  to  greatly 
exceed  the  dynamic  range  limitations  of  the  device  without  causing 
failure  or  damage  to  the  rest  of  the  electrical  system.  Moreover, 
the  device  is  less  likely  to  disturb  the  field  being  measured. 

The  electric  field  generated  by  lightning  strikes  is  generally  much 
lower  than  that  produced  by  EMP,  but  the  duration  is  much  longer. 
Electro  Magnetic  Pulses  typically  have  rise  times  on  the  order  of 
10  ns,  whereas  lightning  strikes  have  rise  times  of  about  6  (is  and 
a  field  strength  of  up  to  1000  V/m  at  a  range  of  1  km. 
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In  a  Mach-Zehnder  interferoneter  without  electrodes,  the  electric 
field  "seen”  by  the  device  will  be  different  froe  the  free  space 
electric  field  because  of  the  pemittivity  of  the  electro-optic 
naterial.  The  device  effectively  acts  as  a  capacitor  in  series 
with  two  free  space  capacitors.  The  dielectric  constant  of  lithiua 
niobate  is  approximately  35.  Thus,  the  field  experienced  by  the 
substrate  is  reduced  by  a  factor  of  about  35. 


The  sensitivity  of  the  device  and  its  dimensions  can  be  determined 
for  a  maximum  electric  field  which  produces  a  peak  deviation  based 
on  linearity  conditions.  The  device  constant  K  is  given  by 


K  = 


(51) 


V,  for  a  push-pull  interferometer  is  given  by 


v;  = 


-^33 


(52) 


for  1  **  1.3vLm,  n.  =  2.15,  r„  =  30.8x10"“, 

L  =  1  mm  and  d  =  1  mm,  V,  »  2000  Volts. 

For  a  phase  deviation  the  voltage  across  the  electro-optic 
material  is 


The  electric  field  inside  the  crystal  is 

„  _  ^LiHbOt 

^LiKbO^  ~  J 


(54) 


The  free  space  electric  field  is  a  factor  of  35  greater  than  that 
in  the  crystal.  Table  VII  illustrates  the  sensitivity  calculations 
based  on  Eqs.  (51)  -  (54).  The  calculations  are  shown  for  a 
crystal  thickness  of  0.5  mm  and  two  different  interaction  lengths 
of  10  and  25  mm.  The  peeUc  deviations  for  two  different  optical 
power  levels  of  1  and  10  mW  are  0.172  and  0.115  radian, 
respectively.  The  electric  field  values  are  that  of  free  space. 

The  25  mm  long  device  is  capable  of  providing  a  sensor  with  a  noise 
floor  of  570  V/m  for  an  optical  power  of  1  mW.  This  device  length 
is  reasonable  for  fabrication  and  will  easily  satisfy  the  frequency 
response  requirements  of  1  GHz.  The  optical  source  power  required 
would  be  2  mW  for  a  lossless  device.  In  practice,  losses  due  to 
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fitwr  coupling  and  scattering  are  about  3  dB.  Hence,  a  acnarce 
of  4  to  5  bN  vould  be  reguired.  Such  laser  scHirces  are 
currently  available  in  the  aarket. 


Tibld  YII 

SEIMTTTVTW  CALmA^»T«iS 


Wavelength  >  1.3|Ui 
Bandwidth  «  IGHz 

CoBpleeentary  detection  using  dual  'i^otodetectors  at  the  output. 
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The  frequency  response  is  lie! ted  priaarily  by  the  transit  tine  of 
the  optical  bean  through  the  device  interaction  length.  The 
crystal's  response  to  a  eicrowave  field  is  instantaneous  on  the 
tine  scale  of  a  microwave  period.  The  microwave  period  must  be 
much  longer  than  this  transit  time  through  the  interaction  region 
for  the  device  to  see  a  static  electric  field.  As  the  microwave 
frequency  increases,  the  optical  beam  sees  a  field  that  is  the 
average  of  the  instantaneous  field  over  the  transit  time.  Were  the 
microwave  period  to  equal  the  transit  time,  the  optical  beam  would 
see  no  net  field,  hence,  no  net  phase  change.  Mathematically,  the 
condition  for  the  upper  frequency  limit  on  the  microwave  signal  may 
be  expressed  as^*'^*: 


2f_T  <  1 


(55) 


where,  f.  »  highest  microwave  frequency 

r  =  transit  time  of  the  optical  beam  through  the 
interaction  region. 


T  is  given  by 


L 

(c/il) 


(56) 


where. 


length  of  the  interaction  region 
velocity  of  light  (3  x  10*  m/sec) 
refractive  index  of  LlNbO,  > 
wavelength) . 


2.15 


Based  on  the  above  equations,  the  highest  frequency  of  operation  is 
calculated  as  follows: 


L  <  70 
L  <  23 
L  <  2.3 


f .  »  1  GHz 
f .  =  3  GHZ 
f.  =  30  GHZ 


The  device  interaction  length  plays  a  major  role  in  determining  the 
fre^ency  response  and  detection  capability.  A  sensor  device 
having  an  interaction  length  of  25  mm  could  operate  from  DC  to 
1  GHz  with  adequate  senstivity  for  a  signal  bandwidth  of  1  GHz. 


^”W.F.  Filter  et.  al.,  "Photonic  measurements  of  microwave 
pulses",  SPIE  pp.  227-232  (1985). 

“J.  Jarzynski  and  R.P.  De  Paula,  "Fiber  optic  electric  field 
sensor  technology",  SPIE  718.  pp.  48-55  (1986). 
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The  extreaely  short  interaction  length  of  2.3  mm  required  to 
achieve  a  sensor  bandwidth  of  30  GHz  would  require  an  optical  power 
of  100  aM  to  obtain  sensitivity  coapareUsle  to  the  3  GHz  device.  An 
alternate  approach  would  be  to  use  the  technique  of  pluuie  reversed 
segaents^^'**  in  which  the  device  bandwidth  could  be  increased 
without  sacrificing  sensitivity. 

We  will  first  describe  the  periodic  phase  reversal  technique  in 
which  the  center  frequency  is  shifted  to  30  GHz.  The  device 
bandwidth  is  typically  Halted  to  less  than  10%,  in  this  exaaple, 
3  GHz.  Figure  30  illustrates  a  Mach-Zehnder  interferoaeter  in 
which  segaents  of  doaain  inverted  regions  are  located  in  the  two 
legs  of  the  interferoaeter.  The  length  of  each  segaent  is 
calculated  based  on  the  center  frequency.  For  a  center  frequency 
of  30  GHz,  the  length  L  is  chosen  to  be  2.3  aa  as  per  the 
calculations  of  Eqns.  (55)  and  (56).  The  phase  reversed  segMnts 
are  separated  by  a  distance  traveled  by  the  optical  signal  in  the 
waveguide  during  one  half  cycle  of  the  aicrowave  signal.  The 
period  of  a  30  GHz  signal  is  33.33  psec,  during  idiich  tiae  the 
optical  signal  travels  4.6  aa  in  a  LiNbO,  waveguide. 

The  device  shown  in  Figure  30  operates  as  follows.  Suppose  the 
aicrowave  signal  is  in  its  positive  cycle  and  a  phase  of  «.  when 
the  optical  signal  enters  the  segaent  A.  As  the  optical  beaa 
travels  through  the  segaents  A^,  and  A«,  it  exi>eriences  equal  and 
opposite  phase  shifts  ±*„  respectively.  The  segaent  spacing  is 
such  that  as  the  optical  signal  enters  segaents  and  B^,  the 
aicrowave  signal  is  in  its  negative  cycle  at  a  phase  of 
Since  the  electric  field  is  now  reversed,  the  segaents  Bt,  and  Bg 
also  experience  ±#^  phase  shifts  respectively.  This  is  because  a 
positive  electric  field  generates  a  positive  phase  shift  in  A,,  and 
a  negative  phase  shift  in  B|,  and  vice  versa.  Hence,  every  segaent 
in  the  lower  and  upper  leg  of  the  interferoaeter  contributes  a 
phase  shift  of  and  respectively.  If  there  are  "n"  segaents 
in  each  leg  of  the  interferoaeter,  the  total  phase  at  the  output  is 

2n*x. 

The  device  shown  in  Figure  30  has  ten  segaents  in  each  leg  of  the 
interferoaeter  for  a  total  length  of  23  aa.  The  segaent  length  of 
2.3  aa  yields  a  center  frequency  of  30  GHz.  The  device  bandwidth 
is  that  of  an  interferoaeter  in  which  the  interaction  length  is 
23  aa,  i.e. ,  3  GHz. 


^Ilod  C.  Alferness  et.  al.,  "Velocity-aatching  techniques  for 
integrated  optic  traveling-wave  switch/aodulators” ,  IEEE  Journal  of 
(2u2mtua  Electronics  OE-20 .  pp.  301-309  (1984). 

**M.  Nazarathy  et.  al.,  "Velocity-aisaatch  coapensation  in 
traveling-wave  aodulators  using  pseudorandoa  switched-electrode 
patterns'*,  J.  Opt.  Soc.  of  Aaer.  A,  pp.  1071-1079  (1987). 
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I  • 


Ulhium  MotMto  Subitrals 


Figure  30.  Periodic  phase  reversal  tecdmique  to  increase 
the  device  operating  frequency  to  30  GBz 
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The  electric  field  seen  by  the  optical  beaa  is  the  average  of  the 
instantaneous  field  over  the  transit  tiae.  If  the  transit  tiae  is 
equal  to  half  the  aicrowave  period,  as  is  the  case  in  this  exaaple, 
the  electric  field  seen  by  the  optical  beaa  is  only  65%  coapared  to 
a  static  electric  field.  This  iaplies  that  the  phase  shift 
experienced  by  the  ten  segaent  interferoaeter  above  is  equivalent 
to  that  of  a  device  having  a  15  aa  long  interaction  region. 

The  periodic  phase  reversed  design  described  above  is  a  narrow  band 
device  that  exhibits  a  bandwidth  of  10%,  i.e.,  3  GHz  for  a  30  GHz 
center  frequency.  A  broadband  device  could  be  designed  in  vhich 
the  segaent  lengths  are  chosen  based  on  a  pseudorandoa  Barker 
code^*.  This  could  yield  highly  sensitive  broadband  devices 
operating  froa  DC  to  40  GHz.  This  design  exhibits  a  non-linear 
phase  response  and  is,  therefore,  applicable  only  in  pulse 
applications. 
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RT.RrTPTr  vjKtn  swusor  calibration 


When  calibrating  and  testing  electric  field  sensors  devised  for 
operation  at  a  few  GHz,  it  will  be  quite  convenient  to  use  snail 
Crawford  TEM  cells.  However,  at  30  GHz,  the  wavelength  in  air  is 
only  1  cm,  and  this  presents  a  problem  in  producing  a  very  compact 
TEM  cell  with  the  desired  uniform  field  characteristics.  The  cell 
would  have  to  be  so  small  that  the  packaged  device  would  not  be 
able  to  fit  into  it.  The  problem  is  with  the  transverse  dimension 
of  the  cell  and  not  with  its  length.  These  cell  dimensions  will 
have  to  be  smaller  than  the  wavelength  in  order  that  there  is  only 
one  mode  within  the  cell  and  that  the  field  is  uniform. 

There  are  two  choices: 

1.  We  can  construct  a  miniature  TEM  cell  cavity  and  test  the 
sensor  in  its  unpackaged  form,  i.e.,  just  the  device  and  its 
attached  fiber  optic  leads. 

2.  Test  the  packaged  device  but  in  free  space.  This  would 
require  the  use  of  an  anechoic  chamber. 

The  first  option  has  the  disadvantage  that  the  package  itself  is 
likely  to  modify  the  response  of  the  sensor.  Hence,  to  test  the 
device  without  its  packaging  may  give  unreliable  results.  The 
second  option  has  the  disadvantage  of  having  to  set  up  a 
transmitting  antenna  that  can  produce  an  accurate  electric  field 
over  a  wide  frequency  range.  It  is  not  inconceivable  that  the 
frequency  calibration  might  be  a  two-step  process.  Over  lower 
frequencies,  a  conventional  TEM  cell  could  be  used,  while  at  high 
frequencies  an  anechoic  chamber  would  be  employed.  This  divide  and 
conquer  approach  would  have  the  advantage  that  the  anechoic  chamber 
would  only  be  used  at  the  highest  frequencies,  and  thus  could  be 
relatively  small. 

Instruments  For  Industry,  IFI,  manufactures  a  Crawford  Cell  TEM 
Test  Chamber,  Model  No.  CCllO,  which  may  be  used  between  DC  and 
1  GHz.  The  first  resonance,  where  the  wavelength  is  equal  to  the 
distance  between  the  top  plate  and  the  septiun,  is  above  1  GHz.  The 
test  chamber  has  inside  dimensions  of  4.7”  X  7.1",  more  than 
adequate  for  the  packaged  sensor.  The  manufacturer  gives  a  maximiim 
recommended  test  object  size  of  2.36”  X  2.36”  X  0.75”  (60  mm  X 
60  mm  X  18  mm) .  The  length  of  the  unpackaged  1  GHz  sensor  is  the 
seune  as  the  maximum  test  object  length.  This  company  could  custom 
build  a  more  compact  TEM  cell  to  take  the  resonance  frequency  to 
about  2  GHz. 

Figure  31  shows  the  three  basic  ways  of  calibrating  the  TEM 
cell.  Method  (a)  uses  an  amplified  wideband  noise  source,  whose 
intensity  is  measured  with  a  spectrum  analyzer  via  the  output  of 
the  electric  field  sensing  system.  Method  (b)  uses  a  tracking 
generator  to  give  a  swept  eunplitude  response  of  the  sensor. 


DAAL02-90-C-0078  79  SRICO 

Optical  Interferometers  for  Sensing  Electromagnetic  Fields 


AmpUnar 


TEMCM 


SOOhmLxMd 


Syatam 


Spwtrum 

Analyzar 


(b)  Tracking  Generator  Method 


(c)  Network  Analyzer  Method 


Figure  31.  Scheaatic  for  calibrating  the  electric  field  sensor 


0AAL02-90-C-0078  80  SRIOO 

Optical  Interforoneters  for  Sensing  Blectronagnetic  Fields 


Method  (c)  uses  a  network  analyzer  to  obtain  the  euiplitude  and 
idsase  response  of  the  sensor. 

Schese  (a)  shows  how  a  wideband  radio  frequency  noise  source  can  be 
used  to  calibrate  the  sensor.  The  output  of  the  device  is 
aaplified  up  to  a  level  that  can  properly  sodulate  the  electric 
field  sensor  via  the  TEM  cell.  The  bandwidth  of  the  device  is 
indicated  by  the  noise  density  profile  (envelope)  as  displayed  on 
the  spectrua  analyzer.  Of  course,  this  scalar  technique  provides 
no  information  about  the  sensor's  phase  response. 

Scheme  (b)  shows  the  preferred  technique  which  uses  a  tracking 
generator  driven  by  the  spectrum  €Ufialyzer.  This  also  is  a  scalar 
measurement  «md  doesn't  yield  any  information  about  the  device's 
phase  response.  The  tracking  generator  is  a  swept  signal  generator 
whose  frequency  is  provided  by  a  local  oscillator  that  is  common  to 
both  it  and  the  spectrum  analyzer.  In  this  way,  at  any  instant  of 
time,  the  frequency  to  which  the  spectrum  analyzer  is  tuned,  is  the 
same  frequency  being  put  out  by  the  tracking  generator.  Thus,  the 
intensity  of  the  received  signal  directly  correlates  with  the 
amplitude  respon.a  of  the  device  under  test. 

One  of  the  most  important  pieces  of  test  instrumentation  in  this 
program  is  the  spectrum  analyzer.  This  will  be  used  for  designing 
and  testing  major  optical  components  in  the  system,  e.g.,  the  laser 
source  and  photodetector,  and  calibrating  the  sensor  in  a  TEM  cell. 
All  aspects  of  the  laser  source  and  photodetector  noise  performance 
can  be  checked  in  the  frequency  domain.  The  tracking  generator 
facility  allows  swept  amplitude  measurements  and,  in  conjunction 
with  another  signal  generator,  gives  us  the  capability  to  perform 
some  limited  two-tone  intermodulation  tests. 

Scheme  (c)  shows  the  more  sophisticated  system  that  cem  yield  the 
phase  response.  The  Network  Analyzer  measures  both  the  amplitude 
and  phase  response.  The  typical  resolution  of  a  microwave  Network 
Analyzer  is  100  kHz.  This  corresponds  to  the  I.F.  bandwidths 
employed  in  these  systems,  and  is  also  about  the  resolution 
bemdwidths  that  will  be  employed  in  the  spectrum  analyzer  systems. 
Such  bandwidths  also  allow  fairly  rapid  sweep  (refresh)  rates  for 
spectrxxm  analyzer  displays  lacking  storage  cap^ilities. 


TKH  CRT  Ji  PBBfPRMftHCB 


The  electric  field  generated  by  this  cheunber  is  given  by: 


(57) 


where  D  is  the  septum  to  top  plate  separation  in  meters.  In  the 
case  of  the  CCllO  TEM  cell,  this  has  a  dimension  of  4.7"/2,  i.e.. 


0AAL02-90-C-0078  81  SRICO 

(^tical  Interferometers  for  Sensing  Electromagnetic  Fields 


2.35"  (60  mm).  To  obtain  a  field  of  1  volt/meter  requires  a  peak 
voltage  of  60  mV,  or  42.4  mV  r.m.s..  This  would  be  produced  by  a 
power  of  36  (iW  or  -14.4  dBm. 


The  graph  shown  in  Figure  32  is  a  plot  of  the  electric  field 
produced  in  the  Crawford  CCllO  TEM  Cell  as  a  function  of  the  drive 
power  to  the  cell,  and  the  phase  modulation  depth  produced  within 
the  cell  for  a  given  drive  power.  The  electric  field,  which  is 
indicated  by  the  full  line,  is  read  off  the  left  hand  ordinate 
scale.  The  induced  phase  deviation,  which  is  indicated  by  the 
broken  line,  is  read  off  the  right  hand  ordinate.  From  Table  VI I, 
Sensitivity  Calculations,  we  see  that  for  a  sensor  with  1  mW 
received  power  at  each  of  the  two  photodetectors  and  a  Nach-Zehnder 
interaction  length  of  25  mm,  the  minimum  detectable  electric  field 
in  a  1  GHz  signal  bandwidth  would  be  570  V/m.  To  produce  this 
field  in  the  TEM  cell  would  require  a  drive  power  of  40.6  dBm, 
i.e.,  11.4  watts.  This  is  positioned  slightly  off  the  right  hand 
edge.  A  drive  power  of  +30  dBm  (1  W)  produces  a  field  of  168  V/m 
and  a  phase  deviation  of  about  1.9  X  10~*  radians.  This  is  about 
10.6  dB  below  the  1  GHz  bandwidth  noise  floor  of  the  system,  and  is 
thus  not  detectable  in  that  bandwidth. 

This  is  not  a  problem  for  the  calibration  procedure  since,  as  we 
have  explained  before,  for  swept  frequency  calibration  purposes  we 
can  employ  spectrum  analyzer  or  network  analyzer  bandwidths  much 
less  than  1  GHz.  If  the  detection  bandwidth  was  100  kHz,  then  the 
signal-to-noise  ratio  and  sensitivity  would  increase  by  a  factor  of 
50  dB.  This  means  that  the  drive  power  requirements  for  a  minimum 
detectable  signal  could  be  reduced  by  50  dBm.  The  TEM  cell  would 
only  require  a  drive  power  of  -9.4  dBm.  In  practice,  we  would 
likely  drive  the  TEM  cell  at  about  +20  dBm  (100  mW)  to  ensure  a 
clean,  low-noise  signal.  If  need  be,  we  could  use  resolution 
bandwidths  as  small  as  1  kHz. 


PYWAMIC.  BAlfgE  TESTS 

In  order  to  exercise  the  device  over  some  dynamic  range,  it  would 
appear  that  we  shall  need  to  use  a  small  RF  power  amplifier  so  that 
powers  of  up  to  about  20  W  (+43  dBm)  may  be  applied.  Based  on 
linearity  considerations,  the  dynamic  range  is  limited  to  a  maximum 
electric  field  of  150,000  V/m.  Clearly,  we  are  not  going  to  be 
able  to  produce  or  sustain  such  a  CW  field  of  this  magnitude,  or 
even  a  pulsed  field.  The  weakest  links  in  the  system  would 
probably  be  the  connectors.  At  high  frequencies,  powers  of  about 
100  W  (+50  dBm)  probably  represent  the  limit.  This  corresponds  to 
electric  fields  in  this  TEM  cell  of  edaout  1.7  kV/m  —  still  very 
small  compared  to  the  150  kV/m  limit  of  the  device  (0.172  radians). 
A  lightning  or  EMP  simulating  test  facility  would  be  required  to 
fully  test  the  device. 
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ELECTRIC  FIELD  &  PEAK  PHASE  MODULATION  DEPTH 
AS  A  FUNCTION  OF  TEM  CELL  DRIVE  POWER 
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Peak  Phase  Modulation  Depth,  radian 


Of  coiirse,  at  low  frequencies  very  high  fields  could  be  produced 
with  close  parallel  plates.  A  voltage  of  150  volts  across  a  1  auD 
gap  will  produce  an  electric  field  of  150  kv/n.  This  should  be  a 
good  way  of  testing  the  dynanic  range  properties  of  the  device  at 
lower  frequencies,  without  concerns  about  the  frequency  response, 
since  the  intermodulation  behavior  should  be  independent  of 
frequency.  A  two-tone  test,  from  two  MHz-type  R.F.  generators 
(frequencies  f^  and  f,)  can  be  used  for  this  test.  R.F.  coils,  or 
transformers,  can  be  used  to  better  impedance-match  the  plate 
intermodulation  test  system  to  the  50-ohm  signal  generators. 
Typically,  f^  »  10  MHz,  and  f,  *  11  MHz,  so  we  would  be  looking  for 
intermodulation  product  frequencies  of  2fx±fa  and  i.e., 
9  MHz,  12  MHz,  31  MHz,  and  32  MHz.  A  100  kHz  I.F.  resolution  would 
easily  discriminate  (separate)  the  9  MHz  and  12  MHz  products  from 
the  10  MHz  and  11  MHz  tones. 

It  should  be  noted  that  in  order  to  the  third  order 
intermodulation  in  a  1  GHz  bandwidth,  we  need  to  apply  an  electric 
field  voltage  consideraby  in  excess  of  150  kV/m.  However,  since  we 
will  be  using  narrow  band  detection,  i.e.,  100  kHz,  we  can  see 
intermodulation  products  50  dB  lower  down  than  would  be  the  case 
for  a  1  GHz  bandwidth,  where  the  dynamic  range,  as  set  by  the 
third-order  intermodulation  is  48.6  dB.  This  is  equivalent  to 
saying  that  the  CNR  is  improved  by  50  dB.  For  small  distortion, 
these  products  are  proportional  to  (4.)’,  so  that  a  50  dB  noise 
floor  reduction  may  be  compensated  by  a  reduction  in  f.  by  a  factor 
of  46.4.  Thus,  the  visibility  of  third  order  intermodulation 
distortion  for  4.  =  0.172  radians  in  a  1  GHz  bandwidth  is  the  same 
as  for  4.  »  3.7  X  10~’  radians  in  a  100  kHz  bandwidth.  At  the  same 
time,  with  its  reduced  deviation,  the  signal  level  has  fallen  by 
16.7  dB.  The  signal  to  distortion  level  in  this  situation  is  now 
about  48.6-16.7+50  »  82  dB.  This  figure  happens  to  be  close  to  the 
limits  on  the  ability  of  most  spectrum  analyzers  to  show  both  the 
signal  and  intermodulation  products  simultaneously  on  the  same 
logarithmic  display.  Of  course,  whether  or  not  both  the  signal  emd 
intermodulation  products  can  still  be  displayed  simultaneously, 
does  not  affect  our  ability  to  measure  each  signal  and  signal 
product. 

We  have  already  seen  that  a  +30  dBm  (1  Watt)  drive  power  to  the  TEM 
cell  produces  a  field  of  168  V/m,  and  a  phase  modulation 
4.  =  1.9  X  10~*  radians.  Thus,  even  the  Crawford  cell  can  manage 
some  limited  intermodulation  tests.  Only  a  test  that  exercises  the 
device  over  w  radian  of  phase  excursion  can  really  test  the 
sensor ^s  linearity.  Its  large  scale  non-linearity  may  not  be 
identical  to  its  small  scale  non-linearity  if  the  cosine-squared 
phase  discriminator  (interference)  characteristic  is  itself  non- 
uniform.  Also,  since  we  would  like  to  see  the  effect  on 
intermodulation  of  bias  drift,  this  can  only  be  done  effectively 
for  large  phase  deviations. 
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APPLICATICBiS 


A  schematic  of  the  sensor  instrument  for  the  remote  measurement  of 
electric  field  or  voltage  is  shown  in  Figure  33.  Several 
applications  for  the  electric  field  sensor  are  listed  below. 


Measurement  of  high  power  microwave  (HPM)  sources  -  The 
sensor  could  be  used  to  detect  leakage  from  industrial- 
grade  radio  frequency  (RF)  heating  and  sealing  systems. 
RF  heating  systems  are  used  in  manufacturing,  food 
processing  and  semiconductor  industries.  Recently,  there 
has  been  an  increasing  public  awareness  of  the  potential 
health  hazard  of  electromagnetic  radiation.  This  issue 
is  presently  undergoing  review  by  the  ERA.  Other 
applications  exist  in  measuring  the  fields  produced  by 
commercial  RF  transmitters  and  antennas  to  monitor  their 
output  and  potential  health  hazard. 

Measurement  of  high  power,  pulsed,  non-repetitive 
microwave  sources  -  Possible  applications  are  detection 
of  high  power  radar  and  electronic  counter  measure  (ECM) 
signals. 

Measurement  of  fields  produced  fay  high  altitude 
electromagnetic  pulse  (HEMP)  simulators  -  Based  on  these 
calculations,  the  insulation  strength  of  high  voltage 
transmission  lines  could  be  designed  to  protect  them  from 
HEMP-induced  flashover.  The  sensor  could  also  be  used  in 
real  time  to  monitor  passage  of  electrical  storms  in  the 
vicinity  of  power  lines. 

Lightning  protection  for  aircrafts  fay  incorporating  the 
sensor  in  their  skins  -  Aircrafts  traditionally  have  used 
metal  skins  which  act  as  a  Faraday  cage  to  protect  the 
electronics  inside  them.  Composite  materials  are 
increasingly  being  used  in  modem  aircrafts  vhich  make 
them  highly  susceptible  to  damage  from  lightning  strikes. 
The  aircraft  can  obtain  advance  warning  of  atmospheric 
conditions  to  avoid  electrical  storms. 

Remote  measurement  of  electric  field  and  voltage  in 
hazardous  and  hostile  environments  -  The  sensor  is 
intrinsically  safe  since  it  carries  no  current  or 
voltage.  Hence,  it  will  not  generate  any  sparks  in 
hazardous  or  explosive  environments.  Since  the  sensor  is 
electrically  isolated,  it  cannot  induce  any  electrical 
breakdown . 
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Near-field  strength  neasurcaent  for  antennas  and  radars 
on  naval  ships  -  The  sensor  could  be  used  to  optimize  the 
transmission  efficiency  and  to  ensure  that  naval 
personnel  are  not  subjected  to  hazardous  levels  of 
electromagnetic  fields. 

Remote  antenna  from  the  shelter  -  The  sensor  could  be 
placed  in  a  forward  position  to  the  receiving  antenna  to 
detect  enemy  counter  measures.  The  information  could 
then  be  transmitted  by  fiber  optic  links  to  a  remote 
shelter  for  automatic  protection  of  sensitive  receivers. 

Medical  and  biological  research  in  the  measurement  of 
electrical  fields  -  The  sensor  could  be  employed  in  the 
test  chamber  to  monitor  and  map  the  fields.  Since  this 
is  a  non-metallic  sensor  it  will  not  disturb  the  fields. 

Measurement  of  electric  field  and  voltage  in  power 
systems  -  The  sensor  could  be  used  in  high  voltage 
environments  which  are  subject  to  a  large  amount  of 
electrical  interference. 
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Figure  33. 


Scbentic  of  a  oonfiguration  for  the  reaote 
aeasureaent  of  electric  field  or  voltage 
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COWCUSIOilS  AMD  RBCOtOmiDATIOIIS 

This  program  has  successfully  demonstrated  the  feasibility  of 
developing  a  nonmetallic  electromagnetic  field  sensor  based  on 
optical  interferometers  in  a  lithium  niobate  electro-optic  crystal . 
This  is  an  all-dielectric  sensor  instrument  that  consists  of 
optical  fibers  for  input  and  output  to  the  electro-optic  crystal. 
Thus,  the  sensor  is  free  from  interference  from  electromagnetic 
fields. 

The  technigue  of  selective  reverse  poling  using  titanium  diffusion 
in  lithium  niobate  was  proposed  and  demonstrated.  This  is  a  novel 
method  to  manufacture  the  sensor  and  SRICO  intends  to  make  an 
application  for  a  United  States  Patent. 

Our  analysis  indicates  that  the  electric  field  sensor  is  feasible, 
and  we  recommend  a  Phase  II  SBIR  program  to  develop  a  prototype 
sensor  instrument.  A  complete  instmimentation  package  including  a 
laser,  a  photodetector  and  signal  processing  electronics  will  be 
developed  during  the  Phase  II  program.  The  sensor  instrument  will 
be  tested  and  delivered  to  the  U.S.  Army  Harry  Diamond 
Laboratories.  We  recommend  the  following  program  to  develop  the 
prototype  sensor  instrument. 

The  objectives  of  this  proposed  Phase  II  effort  are: 

•  to  design  and  develop  an  electrode-less  Mach-Zehnder 
interferometer  for  electric  field  sensing; 

•  to  fabricate  a  prototype  fiber  pigtailed  sensor  device; 

•  to  design  and  develop  a  system  for  testing  and 
calibrating  the  sensor; 

•  to  develop  an  instrumentation  system  incorporating  the 
electric  field  sensor. 

The  overall  program  is  divided  into  two  parts  —  component 
development  and  system  development.  Tasks  1-7,  associated  with 
component  development,  are  designed  to  produce  a  prototype  fiber 
pigtailed  electric  field  sensor.  Tasks  8-10,  associated  with 
system  development,  are  designed  to  produce  the  instrumentation 
system  for  the  electric  field  sensor.  This  Phase  II  SBIR  program 
is  scheduled  over  a  two-year  period. 


Task  1:  Finalize  Electric  Field  Sensor  Specifications 

A  project  kick-off  meeting  will  be  held  with  the  project  manager  to 
finalize  the  performance  requirements  of  the  sensor. 
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TasX  2t _ Finalize  Electric  Field  Sensor  Design 


The  electric  field  sensor  as  per  the  specifications  in  task  1  will 
be  fabricated  in  a  lithium  niobate  substrate  to  operate  at  a  laser 
wavelength  of  1.3  iim.  In  this  task,  the  waveguide  design  will  be 
finalized  for  the  creation  of  the  drawings  to  produce  the 
photomasks.  The  waveguide  will  consist  of  several  test  patterns 
and  the  Mach-Zehnder  interferometer. 

TasX  3; _ Fabricate  Electric  Field  Sensor  for  Tj^hoT-atinry  Testing 

The  sensor  will  be  fabricated  using  the  photomasks  developed  in 
task  2.  Sub-tasks  include: 

Substrate  Preparation  -  The  substrate  purchased  from  the  vendor 
will  be  cleaned  prior  to  domain  inversion  and  waveguide  fabrication 
steps . 

non i n  Inversion  -  Titanivim  metal  strips  will  be  delineated  for 
diffusion  to  produce  selective  domain  inverted  regions  in  the 
crystal . 

Waveguide  Fabrication  -  An  aluminum  metal  mask  will  be  developed  on 
top  of  a  lithium  niobate  crystal  to  act  as  a  barrier  in  the  proton 
exchange  waveguide  process. 

End-Face  Polishing  -  The  end-faces  will  be  polished  for  coupling 
light  in  and  out  of  the  waveguide. 

The  output  of  this  task  will  be  several  waveguide  samples  to  be 
used  for  testing. 

Task  4;  Evaluate  Electric  Field  Sensor 

Light  from  a  1.3  pm  laser  will  be  coupled  to  a  single  mode 
polarization  maintaining  optical  fiber.  The  fiber  output  will  thus 
be  coupled  to  the  optical  waveguide  fabricated  in  the  crystal.  At 
the  output,  a  microscope  objective  will  be  used  to  focus  the  light 
onto  a  photodiode.  The  output  of  this  task  will  be  the 
establishment  of  the  pareuneters  for  the  final  device  design. 

Task  5t  Finalize  Prototype  Electric  Field  Sensor  Design 

Based  on  the  results  of  task  4,  the  prototype  design  will  be 
completed.  In  this  task,  the  waveguide  design  will  be  completed 
for  the  creation  of  the  drawings  to  produce  the  photomasks.  The 
photomask  will  consist  of  the  design  of  single  and  two-ended  Mach- 
Zehnder  interferometer  sensor  devices. 
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TaaK.feS _ Fabricate  Prototype  Electric  Field  Sensor 

The  sensor  will  be  fabricated  using  the  photoaasks  developed  in 
task  5.  The  process  paraaeters  used  to  fabricate  the  sensor  will 
be  Modified  based  on  the  results  of  the  laboratory  device.  Sub¬ 
tasks  include  the  sane  as  those  of  task  3.  Additional  sub-tasks 
are: 

Dice  Crystal  -  The  crystal  sample  will  consist  of  20  or  more  sensor 
devices.  The  crystal  will  be  cut  into  individual  devices  or  dies 
to  facilitate  fiber  attachment  and  packaging. 

Attach  Optical  -  Single  mode  optical  fibers  will  be  attached 
to  the  input  and  output  ports  of  the  sensor  device.  The  input  port 
will  be  attached  to  a  polarization  maintaining  single  mode  optical 
fiber.  The  output  port  will  be  attached  to  stzmdard  communication- 
type  single  mode  fibers. 

The  output  of  this  task  will  be  two  fiber-pigtailed  sensor  devices 
that  will  be  available  for  testing  and  evaluation. 

Tasfc  7; _ Package  Prototype  Electric  Field  Sensor 

The  fiber  pigtailed  sensor  devices  will  be  packaged  in  a  non- 
metallic  enclosure.  The  package  will  be  designed  to  protect  the 
fiber  and  waveguide  interface  from  accidental  damage.  The 
packaging  material  will  be  selected  to  withstand  the  various 
environmental  operating  conditions. 

Task  8;  Design  and  Fabricate  Optoelectronic  Package 

This  task  will  address  the  laser  requirements,  the  electronics  for 
driving  the  laser  and  the  signal  processing  electronics  for 
converting  the  modulated  optical  signal  to  an  RF  signal.  Sub-tasks 
include: 

Laser  Selection  -  The  laser  will  be  selected  based  on  the  optical 
power  required  to  achieve  the  desired  dynamic  range.  Various  laser 
manufacturers  will  be  evaluated  to  obtain  a  low  noise  laser. 
Considerations  will  be  given  to  ease  of  attaching  a  polarization 
maintaining  fiber  to  the  laser  source. 

Laser  Power  Supply  Selection  -  Laser  power  supplies  are  available 
from  several  manufacturers.  The  power  supply  will  be  chosen  based 
on  steUaility  and  its  ability  to  supply  the  required  drive  current. 

Optical  Detector  Selection  -  The  detector  will  be  required  to  have 
a  bandwidth  greater  than  1  GHz.  In  addition,  the  two-ended  Mach- 
Zehnder  interferometer  will  require  a  package  that  has  two 
detectors  in  it. 


DAAL02-90-C-0078  90  SRICO 

C^tuical  Interferometers  for  Sensing  Electromagnetic  Fields 


% 


pjifskftyinff  -  The  optoelectronic  package  will 

be  assenbled  to  contain  the  laser,  the  detector  and  the  associated 
electronics  within  it. 

Fiber  Qptie  mhia  Pahrlfijikion  -  The  fiber  optic  cable  linking  the 
electric  field  sensor  to  the  instrumentation  package  will  contain 
a  single  mode  polarization  maintaining  fiber  and  two  communication 
grade  single  mode  fibers.  If  a  suitable  cable  system  is  not 
available  commercially,  we  will  arrange  to  have  a  special  cable 
made. 

The  output  of  this  task  is  an  instnuentation  package  that  will 
provide  optical  input  to  the  sensor  and  receive  and  process  the 
signal  from  the  electric  field  sensor. 

TaSlL  _ Develop  Calltaration  ■ffar-hnimifts  for  the  Blmctric 

Field  Sfinaor 

A  TEM  cell  will  be  designed  or  a  commercial  unit  modified  and 
interfaced  with  instrumentation  for  evaluating  the  sensor  device. 

Tagfc  IQs _ Evaluate  Prototype  Electric  Field  Eensor 

In  this  task  we  will  evaluate  the  packaged  sensor  device.  The  TEM 
cell  developed  in  task  9  will  be  used  to  evaluate  system 
parameters.  Sub-tasks  include: 

System  Calibration  -  The  device  sensitivity  will  be  measured  to 
obtain  sensor  calibration  curves. 

Moise  Floor  Measnrgment  -  The  minimum  detectable  electric  field 
will  be  measured  with  a  signal  generator  applied  to  the  TEM  cell 
and  a  spectrum  analyzer. 

TnkoiTMvtulation  Distortion _ fnymislc  Ranael  -  By  applying  two 

frequencies,  the  sensor's  dynamic  range  will  be  measured. 

Frequency  Response  Measurement  -  The  TEM  cell  will  be  driven  by  a 
Network  Analyzer  or  a  Tracking  Generator.  Swept  frequency 
techniques  will  be  employed  to  characterize  the  sensor  over  a  1  GHz 
bandwidth.  At  frequencies  below  2  GHz,  a  conventional  TEM  cell 
could  be  used,  while  at  high  frequencies  an  anechoic  chamber  would 
be  required. 

Task  11; _ DgliYfitY_/  Bfiporting 

A  complete  prototype  sensor  unit  will  be  delivered  to  Harry  Diamond 
Laboratories,  at  the  end  of  the  program,  for  testing  and 
evaluation.  Brief  monthly  reports,  a  final  report  and  a  final 
briefing  detailing  the  %K>rk  carried  out  during  this  Phase  II  effort 
will  be  provided. 
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